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Abstract 


Survivability  studies  have  shown  that 

intercontinental  ballistic  missiles  are  vulnerable  to 

thermal  effects.  In  particular,  the  cumulative  thermal 

effect  from  a  multiburst  attack,  and  laser  thermal 

energy  can  seriously  damage  or  destroy  a  missile.  One 

possible  defense  against  the  thermal  threat  is  rotation 

of  the  missile.  The  purpose  of  this  thesis  was  to 

determine  if  rotation’ decreased  the  maximum  skin 

TO 

temperature  of^the  missile  ^increasing  fe-he  missile's 
probability  of  survival. 

The  study  investigated  several  different  scenarios. 

■  v'-.-f  i.v  >• 

The  first  scenario,. was  the  Peacekeeper  Dense  Pack 
missile  system.  The  missile  field  was  subjected  to  a 

walk  attack  of  2  MT  weapons,  with  the  incoming  RV's 

■  \ 

exploding  every  two  seconds.  The(  «-econd/ scenario  was  a 
4-on-l  attack  of  a  missile  launching  system, -Hare 
specifically,  one  missile  was  subjected  to  four  bursts 
located  in  various  positions  surrounding  the  missile. 

The  intent  was  to  determine  if  rotating  a  missile,  even 
when  surrounded  by  thermal  radiation,  would  increase  the 
probability  of  survival.  Finally,  the  missile  is 
attacked  by  a  space-based  laser  with  a  maximum  absolute 


power  of  10  megawatts.  In  all  cases,  the  rotation  rate 
was  limited  to  a  maximum  of  1.6  radians/secwd,  as 
established  by  studies  at  the  Air  Force  Institue  of 
Technology.  Using  computer  programs,  the  maximum  skin 
temperature  was  calculated,  with  the  resultant 
probability  of  damage  determined  using  a  cumulative 
log-normal  distribution  function.  Comparisons  were  made 
between  the  rotating  and  nonrotating  missiles  to 
determine  if  rotation  did  increase  the  probability  of 
survival  for  the  missile  system. 

In  all  scenarios  studied,  rotation  significantly 

.  <  V 

decreased  the  maximum  skin  temperature,  increasing  the  > 
probability  of  sur vival/T or  the  missile.  /The  decrease 
was  most  dramatic  for  the  walk  attack,  where  an  optimum 
rotation  rate  of  .8  radians/secbfld  was  established.  For 
the  4-on-l  attack,  rotation  was  effective,  but  required 
the  maximum  1.6  radians/secb-JKf  rotation  rate  for  best 

.  t. 

results.  Finally,'/  for  the  laser  threat,  rotation  was 
effective  for  the  scenarios  studie dV  with  the  iliaximum 
rotation  rate  providing  the  greatest  amount  of 
protection.  As  a  consequence,  even  at  these  relatively 
low  rotation  rates,  rotation  is  an  effective  defense 
against  the  thermal  threat. 


THERMAL  EFFECTS  ON  A 


ROTATING  MISSILE 

I .  Introduction 

Background 

Thermal  energy  is  electromagnetic  radiation 
travelling  from  a  source  at  the  speed  of  light. 

Normally,  thermal  energy,  such  as  from  the  sun,  produces 
no  harmful  effects.  However,  when  enough  thermal  energy 
is  absorbed  in  a  short  amount  of  time,  it  produces 
heating  and  melting  of  the  absorbing  material.  This 
phenomenon  is  referred  to  as  the  thermal  effect. 

Past  survivability  studies  have  shown  that 
missiles,  because  of  their  thin,  metallic  skin,  are 
vulnerable  to  the  thermal  effect.  This  is  because  the 
missile  surface  absorbs  enough  energy  to  raise  the  skin 
temperature  beyond  the  melting  point,  causing  structural 
failure  and  destruction  of  the  missile.  Two  sources 
capable  of  causing  thermal  damage  to  a  missile  are:  the 
thermal  pulse  from  a  nuclear  explosion,  and  a  laser  beam 
pulse . 

When  a  nuclear  weapon  explodes  near  the  surface  of 
the  earth,  it  creates  a  fireball  emitting  thermal 
radiation.  The  radiation  from  a  single  burst  is  emitted 
as  a  double  pulse,  the  first  pulse  lasting  less  than  a 
tenth  of  a  second,  and  containing  approximately  one 


percent  of  the  energy.  The  second  pulse  may  last 
several  seconds  (up  to  ten  seconds  for  a  one  megaton 
explosion)  and  contains  ninety-nine  percent  of  the  total 
thermal  energy. 

Recent  studies  have  shown  that  when  a  missile  is 
subjected  to  multiple  bursts,  the  thermal  pulses 
overlap,  having  an  additive  or  cumulative  effect  on  the 
missile.  For  these  cumulative  burst  scenarios,  the 
thermal  effect  had  a  significantly  larger  lethal  range 
than  for  multiple  bursts  considered  separately.  This 
lethal  range  was  even  greater  than  the  lethal  range  for 
the  blast  effect.  Thus,  a  missile  is  especially 
vulnerable  to  the  cumulative  thermal  effect  from  a 
multiple  burst  attack. 

If  a  missile  survives  the  nuclear  thermal  threat, 
it  may  still  be  vulnerable  to  the  thermal  effect  from  a 
space-based  laser.  While  the  concept  of  a  laser  in 
space  is  not  new,  the  Strategic  Defense  Initiative  (SDI) 
has  sparked  a  renewed  interest  in  this  type  of  weapon. 
Laser  light  has  several  properties  that  make  it  useful 
for  weapons  application.  First,  the  intensity,  or  power 
per  unit  area,  a  laser  can  deliver  to  a  target  is  very 
large.  Secondly,  laser  energy  is  highly  directional, 
allowing  the  intensities  to  be  confined  to  a  small  area 
and  to  remain  collimated  for  long  ranges.  Thus,  a  laser 
is  capable  of  depositing  a  large  amount  of  energy  in  a 


short  amount  of  time  on  a  relatively  small  area. 

Regardless  of  the  source  of  thermal  energy,  the 
extent  of  damage  to  a  target  depends  on  the  amount  of 
radiation  absorbed  by  a  unit  area  of  missile  surface  in 
a  short  interval  of  time.  For  a  given  surface  material, 
only  a  small  amount  of  the  absorbed  energy  will  be 
dissipated  away  by  conduction,  convection,  or 
re-radiation.  As  a  consequence,  the  absorbed  energy  is 
contained  in  a  shallow  depth  of  the  surface,  resulting 
in  high  temperatures  that  could  damage  the  material. 

The  purpose  of  a  survivability  study  is  to  determine  the 
amount  of  damage  to  the  target  and  express  it  as  a 
probability  of  survival.  For  thermal  effects  on  a 
missile's  skin,  the  probability  of  survival  was 
calculated  using  the  cumulative  log-normal  distribution 
function.  Using  this  technique,  a  probability  of 
survival  was  determined  for  each  missile  as  a  function 
of  the  missile's  maximum  skin  temperature.  A  detailed 
discussion  of  the  distribution  function  is  in  Hall 
(Hall,  1984:  94). 

One  possible  way  to  decrease  the  amount  of  absorbed 
energy,  and  thus  the  skin  temperature  of  the  target,  is 
to  rotate  the  missile.  Rotation  moves  radiated  surfaces 
of  the  missile  out  of  line  of  sight  view  from  the 
source.  The  missile,  in  effect,  acts  as  a  shield 
against  the  thermal  radiation.  Also,  rotation  exposes 


more  surface,  distributing  the  energy  over  a  larger 
area,  keeping  the  skin  temperature  below  the  sure-kill 
level . 

Problem  and  Scope 

The  purpose  of  this  thesis  project  was  to  show  the 
effect  of  missile  rotation  on  the  skin  temperature  of  a 
missile  subjected  to  thermal  radiation.  The  study 
focused  on  the  cumulative  thermal  effects  from  a 
multiburst  attack,  and  the  thermal  energy  deposited  by  a 
hypothetical  space-based  laser  on  a  missile  during 
flight.  The  rotation  rate  was  varied  from  0-1.6  radians 
per  second,  with  a  maximum  skin  temperature  and 
probability  of  survival  calculated  for  each  rotation 
rate.  These  maximum  temperatures  and  probabilities  of 
survival  were  compared  to  the  nonrotating  missiles  to 
determine  if  rotation  increases  survivability.  Also,  an 
optimum  rotation  rate  (one  which  keeps  the  missile  skin 
temperature  the  lowest)  was  determined.  The  intent  was 
to  show  that  since  rotation  distributes  thermal 
radiation  over  larger  areas  of  the  missile's  surface, 
the  thermal  effect  is  reduced,  and  the  missile's 
probability  of  survival  increases. 

For  the  raultiburst  scenarios,  the  burst-target 
missile  system  was  the  Peacekeeper  close-spaced  basing 
(CSB)  system  (see  figure  1).  This  system  was  useful  in 
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modelling  how  a  rotating  missile  would  survive  a  walk 
attack.  Also,  several  4-on-l  attacks  were  studied. 
These  scenarios  involved  an  isolated  missile  subjected 
to  four  bursts,  one  burst  occuring  every  two  seconds. 
The  burst  locations,  relative  to  the  launching  missile, 
were  varied.  The  intent  was  to  determine  if  an  optimum 
rotation  rate  was  scenario  dependent.  The  laser  weapon 
was  a  hypothetical  space-based  design  with  a  maximum 
absolute  power  of  10  megawatts.  While  no  such  system 
exists,  weapons  of  this  type  are  considered  possible. 
Finally,  the  rotation  rate  was  limited  to  1.6  radians 
per  second.  This  maximum  value  was  established  by 
studies  conducted  at  the  Air  Force  Institute  of 
Technology,  Department  of  Aeronautics  (Bandstra,  1985: 


Assumptions  and  Limitations 

The  following  assumptions  and  limitations  were  made 
to  simplify  the  study: 

1.  Existing  computer  programs  were  used  to 
calculate  skin  temperatures  for  this  analysis.  Missile 
flight  characteristics  were  modelled  using  the 
information  in  figure  2.  The  programs  were  modified  to 
account  for  the  missile  rotation,  and  the  missile's 
cylindrical  shape. 
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2.  The  reentry  vehicles  landed  on  their  designated 
targets,  i.e.  no  aiming  error.  This  step  allowed  the 
programs  to  run  on  a  personal  home  computer. 

3.  The  rotation  rate  was  constant.  The  system  was 
assumed  capable  of  reaching  and  holding  the  exact 
rotation  rate. 

4.  The  missile's  cylindrical  shape  was  modelled  as 
a  set  of  evenly  spaced  nodes  of  unit  area.  Three  models 
(8,  16,  and  32  node)  were  evaluated. 

5.  The  laser  weapon  was  a  hypothetical  space-based 
design.  Beam  generation,  power  requirements,  or  other 
beam  propagation  phenomena  were  not  modelled. 

6.  The  laser  beam  spatial  profile  was  assumed  to 
be  flat.  In  otherwords,  the  intensity  in  every  square 
centimeter  of  the  beam  spot  was  assumed  to  be  constant. 
Exact  spatial  profiles  are  hard  to  determine  for  a  high 
power  laser,  and  thus  a  constant  intensity  profile  was 
considered  appropriate  for  this  study  (Bailey,  1985). 

7.  Laser  beam  pulses  were  restricted  to  one  second. 
This  was  considered  an  appropriate  tracking  and  power 
limitation  (Bailey,  1984:41). 

Approach 


The  following  approach  was  used  in  calculating  the 
probability  of  survival  for  a  missile  subjected  to  the 


cumulative  thermal  effects  from  a  multiburst  attack 


Existing  computer  programs  were  used  to  calculate  the 
temperature  on  the  missile's  skin.  These  programs  were 
derived  from  an  energy  balance  over  a  unit  area  of  the 
missile  skin  surface  and  a  simplification  known  as  the 
thin  skin  approximation  (Hall,  1984:12).  The  programs 
were  modified  in  two  ways.  First,  the  temperature  was 
calculated  at  several  locations  or  nodes  around  the 
circumference  of  the  missile,  simulating  the  cylindrical 
shape  of  the  missile.  Second,  the  nodes  were 
repositioned  for  each  time  step,  simulating  the  rotation 
of  the  missile.  In  this  manner,  a  distribution  of 
temperatures  around  the  missile  could  be  estimated.  The 
probability  of  survival  for  the  missile  was  found  by 
using  the  highest  temperature  of  any  node  on  the  missile 
surface . 

The  following  approach  was  used  to  find  the 

probability  of  survival  for  a  missile  subjected  to  a 

laser  pulse.  Considering  the  inherent  limitations  in 

focusing  laser  light,  the  spot  size  varied  from  10  to  40 

cm  in  diameter.  Also,  since  the  laser  pulse  was  limited 

to  one  second  and  the  rotation  rate  was  limited  to  1.6 

radians  per  second,  the  area  of  missile  actually 

irradiated  was  very  small.  Therefore,  the  missile 

surface  was  modelled  as  a  flat  slab,  divided  into  cells 
2 

of  one  cm  area.  Using  an  iterative  process,  the  amount 

9 


of  radiation  absorbed  and  the  temperature  in  each  cell 
was  calculated.  At  the  end  of  the  iteration,  the  cell 
with  the  largest  amount  of  absorbed  energy  and  thus  the 
highest  temperature  was  used  to  calculate  the 
probability  of  survival  for  the  missile. 

Presentation 

Chapter  II  contains  a  derivation  of  the  model  for 
the  cylindrical  shape  of  the  missile,  and  missile 
rotation,  as  well  as  a  discussion  of  the  equations  used 
to  model  laser  energy  deposition.  Then,  how  these 
models  were  incorporated  into  existing  computer  programs 
is  presented.  Chapter  III  summarizes  some  of  the 
parameters  and  conditions  used  in  this  study,  and  why 
certain  parameters  were  chosen.  Finally,  Chapter  IV 
presents  the  results  of  the  study  and  Chapter  V  details 
the  conclusions  and  recommendations. 


II.  Theor 


This  chapter  contains  the  theory  used  to  determine 
the  probability  of  survival  for  a  rotating  missile 
subjected  to  thermal  radiation.  The  first  part  of  this 
chapter  deals  with  the  cumulative  thermal  effect  from 
multiple  nuclear  weapon  explosions.  First,  a  derivation 
of  the  missile's  cylindrical  shape  is  presented, 
followed  by  an  explanation  of  how  the  rotation  of  the 
missile  was  calculated.  The  second  half  of  the  chapter 
examines  laser  radiation.  Specifically,  the  equation 
used  to  calculate  the  amount  of  laser  energy  absorbed  by 
the  missile  surface  is  derived.  Next,  the  laser  spot 
and  missile  surface  are  modelled  to  simulate  rotation. 
Finally,  the  technique  used  to  calculate  maximum 
temperatures  and  probabilities  of  survival  for  the 
missile,  using  these  models,  is  explained. 

Derivation  of  the  Model 
for  the  Missile  *  s  Cylindrical  Shape 

The  missile's  probability  of  survival  for  the 
thermal  threat  from  a  nuclear  weapon  is  calculated  using 
the  maximum  temperature  reached  during  a  burst  scenario. 
The  temperature  is  calculated  over  a  unit  area  of 
missile  skin  surface  using  a  simplification  known  as  the 
thin  skin  approximation.  This  approximation  leads  to 


the  following  differential  equation: 


a 


mi  i 

dt 


incident 

dt 


-  h[T(t) 


T  .  ( t )  ] 
air  '  1 


(2.1) 


which  can  be  solved  using  the  method  of  finite 
differences.  After  making  the  appropriate 
substitutions,  equation  (2.1)  becomes: 
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h  *  t 

[Tj(a-  2maX)  +  h*t 


T  +  a  A  Q  ] 
max  air 


h  *  t 


(2.2) 
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where 


T-  =  temperature  at  the  end  of  the  jth  time  step 

1  (K) 


T,  »  temperature  at  the  beginning  of  the  jth  time 
step  (K) 

a  =  Cppd 

a  »  absorptivity  of  missile  skin 


Q  *  total  thermal  fluence  incident  on  the  missile 
during  the  jth  time  step  ( J/mZ) 


h  «  local  convective  heat  transfer  coeffecient 
(J/m  -s-K) 


T  . 
air 


t 

max 


temperature  of  ambient  air  at  missile  altitute 

(K) 

UU 

time  step  (sec);  t  =.0417*Y*  , 

(Y-yield  in  kilotoni^ 


For  a  more  detailed  derivation  of  this  equation,  refer 


to  Hall  (Hall,  1984:12-18).  Using  equation  (2.2),  the 
temperature  of  the  missile  surface  can  be  found. 

In  past  studies  the  missile's  cylindrical  shape  was 
not  considered,  but  rather  the  missile  surface  was 
modelled  as  a  single  flat  slab.  This  model,  however, 
does  not  adequately  correct  for  the  geometry  between  the 
slab  and  each  burst  point.  For  example,  when  a  missile 
was  subjected  to  more  than  one  burst,  the  thermal 
radiation  from  each  burst  was  assumed  to  strike  the 
slab.  A  correction  factor  was  calculated  for  the 
missile's  flyout  angle,  but  no  correction  was  made  to 
account  for  the  angle  between  the  burst  point  and  the 
slab.  Using  the  cumulative  thermal  energy,  equation 
(2.2)  was  used  to  calculate  the  temperature  of  the 
missile.  In  reality,  the  amount  of  radiation  absorbed 
by  the  slab  also  depends  on  the  geometry  between  the 
slab  and  the  burst.  In  the  extreme  case  of  four  bursts 
completely  surrounding  a  missile,  it  would  be  incorrect 
to  assume  one  spot,  or  slab,  received  all  of  the  thermal 
energy  from  all  four  bursts.  Actually,  the  thermal 
radiation  would  be  distributed  around  the  surface  of  the 
missile,  and  the  resultant  temperature  anywhere  around 
the  missile's  circumference  would  be  much  less  than  the 
temperature  calculated  assuming  one  slab  received  all  of 
the  thermal  radiation.  In  addition,  a  rotating  missile 
would  also  distribute  the  thermal  radiation,  reducing 


the  amount  any  one  point  on  the  surface  absorbed. 


how  much  energy  was  incident  perpendicular  to  the 
surface.  Obviously,  when  the  shape  of  the  missile  is 
considered,  the  surface  facing  the  burst  receives  the 
most  energy,  while  the  surface  on  the  back  receives  no 
energy.  Also,  when  allowed  to  rotate,  the  nodes  could 
be  removed  from  the  side  receiving  radiation  decreasing 
the  amount  of  energy  any  one  node  would  receive. 
Therefore,  to  account  for  the  fact  that  each  node 
receives  a  different  fraction  of  thermal  energy, 
dependent  on  its  position  relative  to  the  burst,  a 
correction  factor  was  needed.  This  rotational 
correction  factor  (RF)  determines  the  fraction  of 
thermal  energy  incident  perpendicular  to  each  node. 
Appendix  A  has  a  complete  derivation  of  the  rotation 
correction  factor.  Now,  with  the  RF  term  included, 
equation  (2.2)  becomes: 


h  *  t 


(n).[V^ - *  h-t.axTair  +  orAQ*RF(n)  ]  (2>3) 


h '  t 


✓  01 3  X  \ 
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where 

n  =  the  node  of  interest,  n=l,2,...8. 

RF(n)=  rotation  correction  factor  for  node  n, 
determines  the  amount  of  thermal  energy 
incident  on  each  node  at  time  step  j. 


Thus,  by  examining  equation  (2.3)  the  temperatures  at 
various  nodes  are  related  to  the  amount  of  energy 
absorbed  by  each  node.  The  value  for  each  rotation 
factor  was  determined  by  its  position  relative  to  the 
burst,  and  rotation  constantly  changed  these  positions. 


Rotation  of  Missile  Nodes 


The  following  definitions  are  needed  to  describe 
the  rotation  of  the  missile: 


tatime  of  missile  skin  exposure  to  thermal 
radiation,  t  is  also  equal  to  the  time 
of  thermal  energy  emission  since  radiation 
travels  at  the  speed  of  light. 


ui  -rotation  rate  of  missile  in  radians  per  second, 
rotation  rate  is  constant,  no  angular 
acceleration 


^t=»time,  in  seconds,  between  succeeding  time  steps 

(  i  .  e  .  it.at..  -  t  .  ) 

J  J  +  l  J 


As  stated  before,  the  cylindrical  shape  of  the 
missile  was  modelled  using  a  set  of  evenly  spaced  nodes 
of  unit  area.  The  position  of  the  nodes  is  determined 
by  their  angular  displacement,  theta  ( 6 ),  and  position 
vector  (r).  Figure  4  shows  an  eight  node  model  with  the 
appropriate  theta  and  r  values. 


kv 


t  of  .59  seconds. 


Theta 

(degrees) 

16.93 

61.93 

106.93 

151.93 

196.93 

241.93 

286.93 

331.93 


Fig.  5  New  Node  Positions 


Thus,  the  new  node  positions  can  be  found  for  each  time 
step,  simulating  the  rotation  of  the  missile. 


Derivation  of  the  Equation 
for  Laser  Energy  Absorption 


Like  the  thermal  radiation  from  a  nuclear 
explosion,  only  a  fraction  of  the  laser  energy  is 
absorbed  by  the  missile  surface.  This  fraction  is  known 
as  the  absorption  coefficient,  a  .  Also,  the  depth  of 
penetration,  or  skin  depth,  5  ,  is  very  small.  Thus,  the 
surface  is  subjected  to  local  heating  and  melting 
( Bai ley , 1 984 : 42 ) .  Appendix  B  contains  sample 
calculations  for  skin  depth,  and  absorption  coefficient 


for  laser  energy. 

Assuming  the  absorption  coefficent  is  known,  the 


rate  of  energy  transfer  into  a  material  surface  is: 


-  a  I  ( J/cm2-s)  (2.5) 

where 
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I  -  radiation  intensity  (J/cm  -s) 
a  ■  absorption  coefficent,  or  absorptivity 

The  total  energy  absorbed  in  a  given  time  interval  will 
be : 


AH  -  t  =  ct  I  At  (J/cm2) 


(2.6) 


If  At  is  replaced  by  a  finite  time  step,  j,  arbitrarily 
chosen  as  .01  seconds,  then  equation  (2.6)  gives  the 
amount  of  energy  deposited  in  each  square  centimeter  of 
the  laser  spot  during  any  time  step.  Knowing  the  amount 
of  laser  energy  deposited  into  the  missile,  equation 
(2.2)  can  be  used  to  find  the  temperature  on  the  missile 
surface . 

Calculating  Laser  Spot  and 
Missile  Surface  Rotation 

For  this  study,  the  laser  spot  size  was  calculated 
using  information  from  Bailey  ( Bai ley ,  1885 : 42-66 ) . 


Using  optimistic  laser  power  and  focusing  capabilities, 

the  spot  size  was  still  quite  small  when  compared  to  the 

missile  surface.  For  this  reason,  the  missile  surface 

was  modelled  as  a  flat  slab.  In  order  to  model  missile 

rotation,  the  missile  surface  was  divided  into  cells, 

2 

one  cm  in  area.  Also,  the  intensity  of  the  laser  was 
assumed  constant  across  the  spot  area.  This  is  a 
reasonable  assumption  for  a  very  high  powered  laser 
(Bailey,  1985).  Thus,  using  equation  (2.6),  the  energy 
deposited  in  each  cell  during  each  time  step  j,  can  be 
calculated.  Since  the  intensity  is  constant,  only  a 
strip  of  cells  were  considered  in  the  calculation. 

Figure  6  shows  a  laser  spot  and  a  strip  of  cells. 


With  no  rotation,  the  spot  would  deposit  the  laser 
energy  on  the  same  cells  for  each  time  interval,  with 
heating  and  possible  melting  of  the  surface.  For 
rotation,  however,  the  spot  would  deposit  energy  on 
different  cells  as  the  missile  surface  moved, 
distributing  the  energy  and  reducing  the  amount  of 
energy  any  one  cell  absorbs. 

The  rotation  of  the  slab  through  the  laser  spot  was 
determined  by  first  calculating  the  linear  velocity  of 
the  surface  of  the  missile  using  the  equation: 

v  a  <D*r  (cm/sec)  (2.7) 

where 

v  *  linear  velocity  of  missile  surface 
(cm/sec) 

a)  rotation  rate  (radians/sec) 
r  a  missile  radius  (cm) 

For  example,  for  a  rotation  rate  of  .8  radians/sec  and 
missile  radius  of  116.84  cm,  the  surface  velocity  is 
approximately  93  cm/sec.  The  distance,  s,  the  slab 
travels  in  any  time  interval  is  found  by: 

s  =  v*At  (cm)  (2.8) 

If  t  is  replaced  by  j,  then  for  each  time  step,  the 
distance  the  surface  moved  in  relation  to  the  spot  is 
known.  The  distance,  s,  determines  which  cells  are 


radiated  by  the  Laser  spot.  The  amount  of  energy  each 
cell  receives  for  each  time  step  is  stored  in  a  one 
dimensional  array.  Obviously,  when  the  laser  spot 
overlaps  a  cell  on  the  surface,  the  energy  in  the  cell 
increases,  causing  higher  temperatures. 

Adapting  Models  t o 
Computer  Programs 

Because  direct  experimentation  is  impossible, 
computer  programs  are  used  extensively  to  model  nuclear 
weapons  effects.  The  groundwork  for  this  study  was 
introduced  in  NE  6.95,  Nuclear  Survivability  of  Systems. 
Also,  Hall  modified  and  improved  the  thermal  model  to 
include  the  cumulative  thermal  effects  from  a  multiburst 
attack  (Hall,  1985:121-134).  This  program,  called 
"Therm",  was  modified  to  account  for  the  missile's 
cylindrical  shape,  missile  rotation,  and  laser  energy 
deposition . 

In  order  to  solve  equation  (2.3),  an  iterative 
method  using  a  finite  time  step  was  used.  Specifically, 
once  the  missile  launch  time,  It,  the  first  burst  time, 
tbl,  and  the  total  number  of  time  steps  were  known,  the 
midpoint  of  each  time  step  was  calculated.  At  these 
times,  missile  characteristics,  ambient  air  temperature, 
and  the  heat  transfer  coefficient  were  calculated  and 


stored  in  one  dimensional  arrays.  In  addition,  using 


equation  (2.4),  new  node  positions  were  found.  Then, 
using  equation  (A. 4),  the  new  rotational  correction 
factors  were  calculated  to  determine  the  amount  of 
radiation  incident  on  each  node.  As  the  temperature  for 
each  node  was  calculated,  a  distribution  of  temperatures 
around  the  circumference  of  the  missile  was  found.  The 
iteration  process  continued,  with  the  maximum 
temperature  achieved  used  to  calculate  the  probability 
of  survival  for  the  missile. 

The  same  computer  program  was  changed  slightly  to 
calculate  the  probability  of  survival  for  a  missile 
subjected  to  a  laser  attack.  Again,  at  the  midpoint  of 
each  time  step,  the  missile  characteristics,  ambient  air 
temperature,  and  heat  transfer  coefficient  were 
calculated  and  stored  in  arrays.  Using  aquation  (2.8), 
the  distance  the  missile  surface  had  moved  was 
calculated.  Representing  the  missile  surface  as  a  one 
dimensional  array,  this  distance  determined  which  cells, 
or  array  locations,  were  radiated  for  that  time  step. 
With  the  appropriate  cells  identified,  the  amount  of 
energy  for  each  cell  was  calculated  using  equation 
(2.6),  and  stored  in  the  array.  At  the  end  of  the 
iteration,  the  cell  with  the  largest  amount  of  energy 
was  used  to  determine  the  maximum  temperature,  and 
probability  of  survival  for  the  missile. 

The  basic  algorithm  for  these  two  models  is  shown 


in  table  I.  In  addition.  Appendix  E.  has  a  complete 
listing  of  the  computer  program  "Therm". 


Il 


TABLE  I 


ALGORITHM  FOR  FINDING  THE  MAXIMUM  SKIN  TEMPERATUE 
FOR  A  ROTATING  MISSILE 


Knowing  It  and  tbl,  calculate  and  store  missile 

velocity,  altitude,  down-range-distance, 

at  t  =  t  +  (j-.5)*T  ,  where 

•io  9i  max* 

j  =  l , 2. . . II , 


Calculate  and  store  T  .  and  h  at  same  times 

air 

Knowing  the  burst  number,  and  the  times 

of  the  burst,  calculate  and  store  the 
amount  of  energy  emitted  by  the  thermal 
source,  and,  if  applicable,  the  height 
of  burst. 

Set  j«l,  T,  a  T  .  at  t  »  t  ,  T~  to  be 

J  j  311*  a  F  2 

any  number  greater  than  T^,  and  the  number 
of  bursts  nb  *  1. 

If  T2  <  T,  and  j  <  =  10  then: 

a.  Tne  current  time  step  is  j  =  j  +  1 

b.  At  time  j*tmax,  determine  if  another  burst 

has  occurred.  If  so,  nb  =  nb  +  1. 

c.  For  each  burst  k  that  occurs: 

1.  Calculate  SR,  theta,  and  thetaprime. 
If  burst  k  has  just  occurred  within 
the  time  step,  missile 
characteristics  must  be 
re-calculated . 

2a.  Calculate  the  total  thermal 
energy  emitted  by  the  burst. 

Calculate  RF(n)  for  each  node, 
and  store  the  amount  of  absorbed 
energy  of  each  node 

2b.  Calculate  the  distance  the  slab  moved 
if  laser  energy  problem,  and  identify 
and  store  the  energy  received  by  each 
cell. 


Calculate  T2  knowing  h;  and  the  total 
energy  absorbed  by  each  node,  or 
each  cell  whichever  is  applicable. 

For  each  node,  or  cell: 

If  T2  <  then  T2  = 

Return  to  condition  in  step  3.  If  either 
test  fails,  go  to  step  4 


Chapter  III.  Conditions  and  Limi ting  Parameters 

In  order  to  calculate  the  probability  of  survival 
for  a  missile  against  any  thermal  threat,  several 
conditions  and  limiting  parameters  were  required.  For 
the  thermal  threat  from  multiple  nuclear  explosions,  the 
conditions  and  parameters  were  the  same  as  those  used  in 
a  similar  study  by  Hall  ( Hall ,  1 985 : 3 1 -34 )  .  A  brief 
description  of  some  of  the  conditions  is  presented 
first,  and  the  remaining  values  listed  in  table  II.  The 
remainder  of  this  chapter  deals  with  the  conditions  and 
parameters  required  to  determine  the  probability  of 
survival  for  a  missile  against  a  laser  threat.  These 
values  were  developed  from  notes  given  by  Dr.  W.  Bailey 
( Bai ley , 1 985 )  ,  and  are  summarized  in  table  III. 

From  a  similar  study  by  Hall,  several  simplifying 
conditions  and  parameters  were  developed  and  used  in 
this  study.  First,  the  RV  aiming  error  was  zero,  thus 
deleting  the  requirement  for  calculating  a  circular 
error  probable.  This  assumption  reduced  computing  time 
significantly,  but  did  not  affect  the  results 
( Ha  1 1 , 1 98 5 : 48 ) .  Second,  the  location  for  calculating 
the  local  heat  transfer  coefficient,  h,  was  assumed  to 
be  the  third  stage  joint  on  a  generic  missile.  This 
position  ( X  =  5.5  m  )  represents  a  point  where  an 
average  amount  of  convective  cooling  occurs  along  the 


missile  skin  (Hall,  1985:34).  Finally,  the  maximum 
number  of  bursts  was  limited  to  four.  While  this  was  a 


rather  arbitrary  limit,  the  cumulative  effects  of  four 
bursts  adequately  demonstrates  the  difference  between  a 
rotating  and  nonrotating  missile.  In  addition, 
increasing  the  number  of  bursts  increases  the  computing 
time,  but  does  not  change  any  other  aspect  of  the 
attack . 

Determining  Minimum  Intensity 
Required  on  Target 

For  this  study,  a  laser  kills  a  target  if  it 

deposits  enough  thermal  energy  per  unit  area  in  a  short 

amount  of  time,  to  raise  the  missile  skin  temperature 

above  the  sure-kill  level  (809  K).  Using  the  thin  skin 
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approximation,  the  minimum  intensity  (W/cm  )  required  to 
insure  a  thermal  kill  can  be  found.  Sample  calculations 
are  in  Appendix  C,  and  they  show  that  approximately  2000 
W/cm  ,  delivered  for  one  second,  will  kill  a  nonrotating 
missile.  Therefore,  for  this  study,  regardless  of  the 
spot  size  considered,  the  minimum  intensity  incident  on 
the  surface  of  the  missile  was  2000  W/cra  . 


TABLE  II 


CONDITIONS  AND  PARAMETERS  FOR  CUMULATIVE  THERMAL  THREAT 


System : 

Close-Spaced  Basing,  of  Peacekeeper  Missiles 
(see  figure  1) 

Threat  Conditions; 

Walk  attack  starting  at  silo  #1  and  continuing 
every  2  seconds  on  successive  silos 
Weapon  Yield:  2  MT 
Height  of  Burst:  0  m 

For  surface  bursts,  tf  =  .18 
(Glasstone  and  Dolan,  1977:319) 

Missile  Conditions : 

Missile  velocity,  altitude,  down-range-distance, 
and  flight  path  angle  shown  in  figure  2  as 
a  function  of  time 
Skin  material:  Aluminum 
K  =  .0001  n>  /s 
p  a  27 00  kg/mJ 
«  =  .50 

C  a  900  J/kg-K 
I  p  a  619  K 
I*®  =  809  K 

Ski nthickness :  d  =  .001  m 
Rotation  Rate:  0  -  1.6  radians/s 


Probability  Conditions: 

RV  aiming  error:  none 

Probability  of  damage  based  on  intensity  and 

calculations  using  the  cumulative  log-normal 
distribution  function 

Pd(Iss)  -  -98 

'd'O  -  -02 


Parameters : 

Maximum  number  of  bursts  considered:  maxb  =  4 

Maximum  number  of  time  steps  needed:  11 

Heat  transfer  coefficient  calculated  at  x  a  5.5  m 


TABLE  III 


CONDITIONS  AND  PARAMETERS  USED  IN  LASER  THREAT  STUDY 


System: 

Peacekeeper  Intercontinental  Ballistic  Missile 

Threat  Conditions: 

Laser:  Space  based,  continuous  power, 
chemical  laser 
Maximum  Power  =>  10  MW 
Aperture  =>  10  m 
Altitude  =*  400-800  km 
Spot  size  -10-40  cm 
Minimum  Intensity  required  on  target: 

2000  W/cra  incident  normal  to  surface 
Maximum  pulse  time  -  1  second 

Missile  Conditions: 

Velocity  =  5  km/s 

Altitude  -  100  km 

Skin  material:  Aluminum,  properties  in  table  II 

Skin  thickness:  d  »  .001  m 

Rotation  Rate:  0  -  1.6  radians/s 

Probability  Conditions: 

Same  as  table  II 

Parameters : 

Maximum  number  of  pulses  on  target:  1 


Determining  Laser 
Threat  Conditions 


As  stated  before,  the  laser  weapon  was  assumed  to 

have  a  maximum  absolute  power  of  10  MWs,  with  an 

aperture  of  10  m.  While  these  are  optimistic 

parameters,  the  possible  deployment  of  such  a  weapons  in 

the  future  is  considered  feasible  (Bailey,  1985).  In 

order  to  calculate  the  probability  of  survival  for  a 

missile  subjected  to  a  laser  pulse,  a  specific 

laser-missile  scenario  was  developed  from  notes  given  by 

Dr.  W.  Bailey  (Bailey,  1985:  42-62). 

The  laser's  initial  intensity  is  given  simply  by 

the  equation:  Intensity  =  Power/Area.  For  this 
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particular  laser,  I  *  1.27*10  W/cm  .  To  be  a  lethal 
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beam,  it  must  be  focused  to  a  small  spot  on  the  target. 
The  spot  radius  for  a  given  distance,  Z,  is  determined 
by : 


W(z) 


.6*A*Z 


(3.1) 


where 

A  =  laser  wavelength  (m) 


Z  =  distance  from  laser  to  target 
=  laser  aperture  radius  (5m) 


Table  IV  shows  several  values  of  Z,  and  the  appropriate 
spot  radii. 
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TABLE  IV 

SPOT  RADII  FOR  VARIOUS  DISTANCES  FROM  LASER  TO  TARGET 


Distance 

(km) 

250 

420 

625 

800 

2000 


Spot  Radius 
(cm) 

10 

20 

30 

40 

100 


Assuming  no  beam  jitter,  the  intensity  on  target  is 
defined  as: 


W. 


'I 

Ifc  «  IQ  u(z) 2  [absorption  +  scatter]  (3.2) 


where 


2 

I  *  intensity  on  target  (W/cro  ) 


2 

I  =  initial  intensity  of  laser  (W/cm  ) 


Wx  *  laser  aperture  radius  (5m) 


W(z)  a  spot  radius  (cm) 


According  to  Bailey,  for  the  altitudes  considered  in 
this  scenario,  absorption  and  scatter  of  the  beam  is 
negligible  (Bailey,  1985:9-11).  Therefore,  the 
intensity  on  target  is  simply: 
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v%s 


(W/ cm  ) 


(3.3) 


I 


t 


!o  W(  z ) 2 


The  maximum  spot  size,  and  thus  the  highest 
altitude  the  laser  could  be  stationed  and  still  deliver 
a  lethal  intensity  can  be  found  using  equation  (3.3)  and 
solving  for  W(z).  Rearranging  the  equation: 


W(z) 


W 

r  I  1  1.5 

[  O -  ] 

I 


(m) 


(3.4) 


where,  for  this  specific  case, 

I  =  1.27  *  105  W/m2 
o 

It  =  2.00  *  107  W/m2 
Wx  =  5m 

W(z)  *  maximum  spot  size 


W(z)  is  approximately  40  cm.  This  means  that  a  10  MW 

laser,  attacking  targets  approximately  800  km  uway,  will 

produce  a  spot  with  a  40  cm  radius,  and  deposit 
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approximately  2000  J/cm  on  the  target  in  one  second. 

This  is  sufficient  to  kill  the  target.  For  closer 

distances,  the  power  required  to  produce  the  desired 

intensity  is  decreased  as  the  spot  size  decreases. 

However,  for  this  study,  the  intensity  on  target  was 

2 

assumed  to  remain  constant  at  2000  W/cm  ,  regardless  of 
the  spot  size  being  considered. 


Chapter  IV  Results  and  Discussion 
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The  main  objective  of  this  thesis  was  to  determine 
whether  rotation  increased  the  probability  of  survival 
for  a  missile  subjected  to  the  thermal  threat.  The 
results  from  the  rotating  case  were  compared  to  the 
nonrotating  case  for  both  sources  of  thermal  energy: 
nuclear  weapons  explosions,  and  laser  pulses.  The 
comparison  showed  that  rotating  missiles  experience  a 
significant  decrease  in  the  maximum  skin  temperature 
regardless  of  the  source  of  the  thermal  energy.  Thus, 
rotation  is  an  effective  defense  against  thermal 
radiation. 

The  results  for  a  missile  subjected  to  the 
cumulative  thermal  effect  from  multiple  bursts  of  a  walk 
attack  are  presented  first.  In  all  cases  the  time  step 
was  tmax/2.  As  a  result,  the  maximum  temperature  change 
between  time  steps  was  less  than  seven  degrees.  Thus, 
accurate  temperature  rises  were  calculated  while  keeping 
the  computing  time  to  a  minimum.  To  begin,  a  comparison 
between  the  8,  16,  and  32  node  models  was  made  to 
determine  which  model  best  simulated  the  missile's 
cylindrical  shape.  Next,  the  effect  of  rotation  on  the 
maximum  temperature,  and  probability  of  survival  of  a 
missile  are  given.  These  results  show  that  an  optimum 
rotation  rate  exists.  Next,  the  scenario  is  extended  to 
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show  the  sure-kill  and  sure-safe  regions  for  the  entire 
missile  field  when  every  missile  is  rotating  at  the 
optimum  rotation  rate.  In  addition,  the  results  for  a 
missile  subjected  to  a  4-on-l  attack  are  presented. 

Under  this  type  of  attack,  a  faster  rotation  rate  is 
desirable.  Finally,  a  comparison  is  made  betweeen  the 
blast  effect,  and  the  thermal  effect  when  rotation  is 
considered.  These  results  show  that  rotation  decreases 
the  sure-kill  region  significantly,  reducing  the 
devastating  effect  of  cumulative  thermal  radiation. 

In  a  similar  manner,  the  effect  of  rotation  on 
maximum  temperature  and  probability  of  survival  are 
presented  for  the  laser  energy  threat.  This  comparison 
shows  that  for  the  scenarios  investigated,  rotation  can 
significantly  increase  the  probability  of  survival  for  a 
missile.  Finally,  the  effect  of  spot  size  on  rotation 
is  examined,  emphasizing  that  large  spot  sizes,  when 
coupled  with  sufficiently  high  energy  densities,  are 
difficult  to  defend  against. 


Comparison  between  8.,  1 6 , 
and  32  Node  Models 


Figure  7  shows  the  probability  of  survival  versus 
rotation  rate  for  missile  #41  using  an  8,  16,  and  32 
node  model.  The  16,  and  32  node  models  agree  over  the 
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entire  range  of  rotation  rates,  but  the  8  node  model 
oscillates  in  the  lower  rotation  rates.  Computing  times 
did  increase,  but  not  significantly.  Therefore,  to 
insure  accuracy  while  keeping  the  computing  time  to  a 
minimum,  the  16  node  model  was  chosen  for  all  subsequent 
results . 

Effect  of  Rotation  on  Temperature 
and  Probability  o f  Survival  for 
Multiburst  Case 

Figure  8  shows  the  maximum  temperature  of  any  node 
versus  rotation  rate  for  missile  #41  subjected  to  four 
bursts.  The  missile  was  launched  at  the  same  time  that 
silo  #1  was  hit,  and  fireball  rise  was  considered. 
Tabulated  data  for  the  curve  is  in  Appendix  D. 

The  effect  of  rotation  on  maximum  temperature  is 
clearly  illustrated  in  figure  8.  As  the  rotation  rate 
increases,  the  temperature  decreases  to  a  minimum,  and 
then  remains  relatively  constant,  oscillating 
approximately  20  degrees  above  the  absolute  minimum. 

The  minimum  rotation  rate  correspondes  to  a  period  of 
rotation  of  approximately  eight  seconds.  For  faster 
rotation  rates  the  nodes  return  to  their  beginning 
positions  too  quickly,  causing  higher  temperatures. 
Another  minimum  is  seen  at  approximately  twice  the 
period  (1.6  radians/sec).  The  drastic  temperature 
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decrease  occurs  because  the  rotation  constantly  moves 
the  nodes,  shielding  them  from  radiation,  while  exposing 
new  nodes  (missile  surface)  to  the  radiation.  This 
prevents  any  one  node,  or  point  on  the  surface  of  the 
missile  from  receiving  too  much  radiation.  Also,  when 
shielded,  convective  cooling  can  be  more  effective, 
allowing  the  nodes  to  reach  lower  temperatures  before 
being  reradiated  on  their  return  trip. 

Figure  9  shows  the  probability  of  survival  versus 
rotation  rate  for  the  same  missile.  This  figure  shows 
how  increasing  the  rotation  rate  increases  the 
probability  of  survival.  However,  the  probability  of 
survival  is  one  at  a  much  lower  rotation  rate  (.4 
radians/sec)  than  the  rotation  rate  where  the  minimum 
temperature  occurs  (.8  radians/sec).  This  happens 
because  the  probability  of  survival  is  calculated  using 
a  cumulative  log-normal  function.  Thus,  once  the 
sure-safe  temperature  is  reached,  the  probability  of 
survival  is  essentially  one.  Therefore,  the  optimum 
rotation  rate  was  defined  to  be  the  rotation  rate  that 
produced  the  lowest  missile  skin  temperature.  For  all 
missiles  in  the  close-spaced  base  system,  the  optimum 
rotation  rate  was  .8  radians  per  second. 
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Comparing  the  Rotating 
and  Nonrotating  Cases 

Table  IV  shows  a  comparison  of  maximum  temperature 
and  probability  of  survival  for  several  missiles  at  no 
rotation,  and  at  the  optimum  rotation  rate. 


TABLE  IV 


COMPARISON  OF  MAXIMUM  TEMPERATURES  AND  PROBABILITY 
OF  SURVIVAL  FOR  NO  ROTATION  AND  OPTIMUM  ROTATION 


Missile  launch  time: 

0  sec 

Time  of  first  burst: 

0  sec 

Time  between  bursts: 

2  sec 

1  cell  CEP 

Optimum  rotation  rate:  .8 

radians/ sec 

Missile# 

No  Rotation 

MaxT  (K)  Ps 

Optimum  Rotation 
MaxT  (K)  Ps 

28 

1638 

0.00 

843 

0.00 

29 

1390 

0.00 

735 

0.28 

30 

1383 

0.00 

731 

0.31 

31 

1169 

0.00 

645 

0.92 

32 

1178 

0.00 

645 

0.92 

33 

1160 

0.00 

640 

0.93 

34 

1015 

0.00 

578 

0.98 

35 

1011 

0.00 

577 

0.98 

36 

885 

0.00 

527 

0.99 

37 

889 

0.00 

527 

0.99 

38 

880 

0.00 

524 

0.99 

39 

789 

0.05 

486 

1.00 

40 

787 

0.05 

486 

1.00 

41 

708 

0.49 

454 

1.00 

42 

680 

0.50 

440 

1.00 

43 

670 

0.54 

436 

1.00 

44 

645 

0.92 

430 

1.00 

45 

640 

0.93 

430 

1.00 
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For  this  thesis,  the  sure-safe  limit  is  defined  by  a 
probability  of  damage  of  .02,  and  the  sure-kill  limit  by 
the  probability  of  damage  of  .98.  Any  values  lower  than 
.02  are  rounded  to  0,  and  any  values  higher  than  .98  are 
rounded  to  1.  Any  missile  not  listed  on  table  IV  has  a 
probability  of  survival  of  0  or  1,  depending  on  the 
missile's  position. 

Another  way  of  presenting  the  values  given  in  table 
IV  is  by  using  figure  1  to  illustrate  the  sure  safe  and 
sure-kill  regions.  Figure  10  shows  the  two  regions  for 
the  values  in  table  IV.  Comparing  these  results  shows 
that  rotation  significantly  decreases  the  sure-kill 
region  for  thermal  effects. 

Comparing  Cumulative  Thermal  Effects 
to  Blast  Effects 

According  to  a  recent  study,  the  cumulative  thermal 
effect  from  a  multiburst  attack  is  more  lethal  than  the 
( noncumulati ve )  blast  effect  (Hall,  1985:52).  Table  V 
shows  the  probability  of  survival  for  noncumulati ve 


blast  effects . 


TABLE  V 

RESULTS  FOR  NONCUMULATI VE  BLAST  EFFECTS 


Missile  launch  time:  0  sec 
Time  of  first  burst:  0  sec 
Time  between  bursts:  2  sec 


Missile  # 

Prob.  of  Survival 

23 

0.000 

24 

0.072 

25 

0.096 

26 

0.590 

27 

0.626 

28 

0.750 

29 

1.000 

Comparing  these  results  to  table  IV  shows  that  the  blast 
effect  is  overwhelmed  by  the  cumulative  thermal  effect 
for  nonrotating  missiles.  However,  for  missiles 
rotating  at  the  optimum  rotation  rate,  the  sure-safe  and 
sure-kill  regions  are  decreased  significantly.  Figure 


11  shows  these  regions  for  rotating  and  nonrotating 
missiles  compared  to  the  blast  effect.  If  cumulative 
blast  effects  were  considered,  the  gap  between  thermal 
and  blast  sure-kill  and  sure-safe  regions  would  probably 
be  decreased  even  more.  These  results  emphasize  how 
rotation  reduces  the  dominance  of  the  cumulative  thermal 


effect  from  nuclear  weapons  explosions, 


Probability  o f  Survival 
and  Maximum  Temperature  for 
4-on- 1  Attack 


■>: 


The  previous  results  concerned  the  thermal  threat 
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Fig.  11.  Location  of  Sure-kill  and  Sure-Safe  Regions 


from  a  walk  attack,  where  the  bursts  progress  in  an 
orderly  fashion  (for  example,  south  to  north).  In  the 
4-on-l  attack,  the  missile  was  subjected  to  four  bursts. 
The  burst  locations,  relative  to  the  launching  missile, 
were  such  that  the  missile  was  radiated  from  all  sides. 
The  intent  was  to  determine  if  the  optimum  rotation  rate 
was  scenario  dependent.  Each  burst  occurred  every  two 
seconds,  and  fell  in  a  counterclockwise  manner  about  the 
launching  missile. 

Figure  12  shows  the  probability  of  survival  vs. 
rotation  rate  for  a  missile  under  a  4-on-l  attack.  The 
probability  of  survival  initially  goes  down,  reaching  a 
minimum  at  approximately  .7  radians/sec.  This  occurs 
because  the  missile  rotation  rate,  and  the  burst 
explosion  rate  are  approximately  the  same.  In 
otherwords,  the  missile  was  rotating  into  the  thermal 
pulse  as  each  weapon  exploded.  However,  as  the  rotation 
rate  increased,  the  probability  of  survival  also 
increased.  Thus,  at  the  higher  rotation  rates,  the 
thermal  energy  was  effectively  distributed  on  the 
missile  surface,  decreasing  the  maximum  temperature,  and 
increasing  the  probability  of  survival.  Several  other 
4-on-l  scenarios  were  examined,  and  in  all  cases,  the 
maximum  rotation  rate  was  the  most  effective. 


Effect  of  Rotation  on 
Temperature  for  the  Laser  Threat 

Figure  13  shows  the  maximum  temperature  versus 

rotation  rate  for  the  missile  surface  subjected  to  a 

laser  pulse.  The  laser  spot  radius  on  target  was  10  cm, 

2 

with  an  intensity  of  2000  J/cm  -s. 

The  effect  of  rotation  on  maximum  temperature  is 
clearly  shown  for  this  scenario.  As  the  rotation  rate 
increases,  the  maximum  temperature  decreases  to  a 
minimum.  The  high  temperatures  at  the  low  rotation 
rates  are  a  result  of  the  laser  spot  overlapping  on  the 
missile  surface  during  the  iteration  process.  As  the 
rotation  rate  increases,  the  spot  overlaps  less  until 
for  each  time  step,  the  laser  spot  moves  an  entire  spot 
width,  and  no  overlapping  occurs.  As  the  rotation  rate 
increases,  and  the  temperature  decreases,  the 
probability  of  survival  increases  accordingly.  Thus, 
for  this  limited  scenario,  rotation  significantly 
increases  the  probability  of  survival.  However,  spot 
size  drastically  changes  the  rotation  rate  at  which  the 
minimum  temperatures  are  reached. 
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Figure  14  shows  the  probabilty  of  survival  versus 

rotation  rate  for  a  missile  subjected  to  a  one  second 
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laser  pulse  with  an  intensity  of  2000  J/cm  -s.  The  four 
curves  are  for  the  four  spot  sizes  considered:  10,  20, 
30,  and  40  cm  radius.  The  most  obvious  characteristic 
of  the  illustration  is  that  as  the  spot  size  increases, 
the  rotation  rate  required  to  cause  an  increase  in 
survival  also  increases.  This  seems  logical  since  a 
larger  spot  would  require  a  higher  rotation  rate  to 
reduce  the  spot  overlap.  For  rotation  rates  restricted 
to  1.6  radians  per  second,  a  large  spot  could  not  be 
moved  quickly  enough  to  reduce  the  missile  skin 
temperature.  In  that  event,  rotation  would  have  little 
effect  on  the  probability  of  survival  for  the  missile. 


Chapter  V.  Conclusions  and  Recommendations 

Conclusions 

From  the  results  given  in  Chapter  IV,  the  following 
conclusions  are  made: 

1.  The  maximum  skin  temperature  for  a  missile 
exposed  to  thermal  radiation  can  be  decreased  if  the 
missile  rotates.  The  most  effective  rotation  rate 
depends  on  the  scenario.  However,  in  general,  the 
faster  the  rotation  rate,  the  greater  decrease  in  skin 
temperature.  Therefore,  rotation  is  an  effective 
defense  against  the  thermal  threat. 

2.  In  a  scenario  such  as  a  walk  attack,  where  the 
attack  progresses  from  one  quadrant  (i.e.  south  to 
north),  rotation  can  have  dramatic  effects.  In 
particular,  a  relatively  low  optimum  rotation  rate  of  .8 
radians/sec  was  observed.  Using  this  optimum  rotation 
rate,  the  sure-kill  region  for  the  cumulative  thermal 
threat  can  be  significantly  decreased.  Specifically, 
for  the  walk  attack  scenario,  10  more  missiles  survive, 
and  the  sure-kill  region  is  reduced  by  approximately 
2000  m. 

3.  For  the  4-on-l  attack,  where  four  bursts 
surround  a  launching  missile,  rotation  still 
significantly  decreases  the  maximum  skin  temperature. 
However,  the  most  effective  rotation  rate  for  this  type 


of  attack  was  the  maximum  rate  of  1.6  radians/sec. 

4.  As  stated  before,  the  cumulative  thermal  effect 
is  more  lethal  than  the  noncumulati ve  blast  effect. 
However,  when  rotation  is  considered,  the  dominance  of 
the  thermal  threat  is  decreased  to  a  point  where  the 
sure-kill  regions  for  the  blast  effect  and  the  thermal 
effect  are  approximately  the  same. 

5.  Rotation  is  an  effective  defense  against  the 
laser  energy  threat.  For  the  scenarios  considered,  the 
1.6  radian/sec  rotation  rate  was  adequate  to  keep  the 
missile  skin  maximum  temperature  below  the  sure-kill 
level.  However,  rotation  as  a  defense  against  the  laser 
threat  does  have  certain  limitations.  A  sufficiently 
powerful  laser,  with  a  large  spot  size  could  kill  a 
missile  before  rotation  could  remove  the  missile  surface 
from  the  spot.  Therefore,  for  rotation  to  be  effective, 
the  threat  against  the  missile  must  be  known. 

Recommendations 

Based  on  the  assumptions  presented  in  Chapter  I,  as 
well  as  the  observations  made  during  the  study,  the 
following  recommendations  are  made: 

1.  Since  the  effectiveness  of  rotation  is  scenario 
dependent,  a  more  thorough  investigation  should  be  made 
into  different  scenario  types. 

2.  A  study  should  be  made  using  the  more  accurate 


cylindrical  geometry  between  the  missile  surface  and 
laser  spot.  This  is  especially  true  as  the  laser  spot 
size  increases  to  a  point  where  the  entire  diameter  of 
the  missile  is  covered  by  the  spot. 

3.  As  recommended  in  a  previous  study  by  Hall,  the 
synergistic  effects  of  thermal  heat  and  blast  should  be 
studied  (Hall,  1985:66).  This  study  did  not  investigate 
the  mechanical  stresses  imposed  on  a  rotating  body. 
However,  the  increased  load  factor  caused  by  rotation, 
coupled  with  the  heat  and  blast  from  a  nuclear 
explosion,  may  exceed  the  structural  limits  of  the 


missile. 


Appendix  A.  Calculating  the  Rotation  Correct  ion  Factor 

This  appendix  outlines  the  method  used  to  calculate 
the  rotation  correction  factor  (RF).  This  correction 
factor  accounts  for  the  cylindrical  shape  of  the 
missile.  A  correction  factor  was  needed  to  determine 
the  amount  of  thermal  radiation  that  falls  incident 
perpendicular  to  the  missile  skin's  surface.  An  8  node 
model  will  be  used  for  this  illustration. 

The  Rotation  Correction  Factor 

The  RF  was  found  by  considering  the  geometry 
between  the  burst  point  and  the  8  nodes  that  make  up  the 
missile's  cylindrical  shape.  For  this  study,  the  RV '  s 
were  assumed  to  land  directly  on  their  targets,  with  no 
aiming  error.  Thus,  the  burst  locations  are  known 
exactly.  The  next  step  was  to  find  the  exact  location 
of  the  8  nodes. 

Using  a  standard  polar  coordinate  system,  the 
location  of  the  8  nodes  was  known  for  each  time  step  by 
their  position  vector,  r,  and  angular  displacement, 
theta  (0).  Since  the  missile  radius  is  very  small 
compared  to  the  distance  between  the  burst  and  missile, 
only  the  angular  displacement  was  used.  The  next  step 
was  to  find  the  relative  position  where  the  thermal 
radiation  is  perpendicular  to  the  missile  surface.  This 
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Figure  A.l  Sample  Burst-Missile  Encounter 
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thetaprime,  the  rotation  correction  factor  for  each  node 
can  easily  be  determined.  First,  recognizing  that  the 
cosine  of  the  difference  between  thetaprime  and  the 
angle  theta  for  each  node  is  the  measure  of  the  fraction 
of  the  thermal  radiation  perpendicular  to  that  node,  the 
RF  for  any  node  can  be  found  using  the  following 
equation  : 

RF(n)  =  Cos (The ta prime-theta ( n  )  )  (A. 4) 

where 

thetaprime  =  location  on  missile  surface  where 
radiation  is  perpendicular 

theta(n)  =  location  of  node  (n  =  1,2, ...8) 

RF(n)  =  correction  factor  for  node,  n 

For  each  time  step,  theta(n)  changes,  which  in  turn 
changes  the  value  of  the  correction  factor  RF(n).  When 
the  difference  between  the  two  angles  is  greater  than  90 
degrees,  the  cosine  becomes  negative.  This  means  the 
radiation  is  shielded  from  the  node,  and  therefore,  the 
RF(n)  is  set  to  zero.  This  of  course  assumes  perfect 
shielding,  and  does  not  account  for  scattering  of  the 
radiation.  However,  radiation  scattering  was  not 
modelled  in  this  study. 
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Appendix  B.  Skin  Depth .  and  Absorption  Coef f icient 

This  appendix  shows  sample  calculations  of  skin 
depth  and  absorption  coefficient  for  laser  radiation  on 
aluminum.  The  laser  is  assumed  to  be  a  chemical  laser 
with  an  operating  wavelength  of  4.0  *  10"^  m.  Further 
examples  can  be  found  in  Bailey  (Bailey,  1985:  42-44). 


Skin  Depth 

Whenever  electromagnetic  radiation  strikes  a  good 
conductor,  like  aluminum,  it  can  be  absorbed  or 
reflected.  In  either  case,  the  electromagnetic  wave 
must  penetrate  the  material  in  order  to  interact  with 
the  material's  electrons.  The  depth  of  penetration  can 
be  found  from: 


<  *c*  A  .5 

rr  rr  J 


(B.l) 


where 

(t  =  3*10  mho/m 
K  m  4*10-6  m 
t  =  9*10~12  C2/ J-m 
c  »  3*10^  m/sec 


Substituting  into  equation  (B.l),  skin  depth  is 

—  8 

approximately  1.1*10”  m.  Thus,  the  heating  is  confined 


to  the  surface  of  the  aluminum  skin. 


Absorption  C o e f fici e n t 


The  fraction  of  energy  actually  absorbed  by  a 
material  is  known  as  the  absorption  coefficient.  For 
good  conductors,  it  is  found  using  the  equation: 


a 


4*^  8 
[  A 


] 


(  B  .  2  ) 


Substituting  into  this  equation,  alpha  is  approximately 
.03.  However,  this  calculation  is  based  on  an  ideal 
material  surface.  The  actual  value  is  expected  to  be 
higher,  because  the  missile  surface  has  impurities, 
oxidation,  and  defects  which  would  increase  the  energy 
coupling  to  the  surface.  Ready  suggests  a  value  of  .1 
( Ready , 1971 : 49 ) .  Therefore,  for  this  study,  the 
absorption  coefficient  is  .1  for  all  laser  calculations 


Appendix  C.  Calculation  of  Minimum  Energy  Density 

The  minimum  energy  density  required  to  achieve  a 
thermal  kill  is  calculated  using  the  thin  skin 
approximation.  Specifically,  the  amount  of  absorbed 
thermal  energy  required  to  raise  the  missile  skin 
temperature  above  the  sure-kill  level  of  809  degrees  in 
one  second,  is  determined.  Also,  this  calculation  was 
checked  by  two  other  methods.  However,  only  the  results 
of  these  methods  are  presented,  along  with  appropriate 
references . 

Calculating  Minimum 
Energy  Density 

As  explained  in  chapter  II,  the  missile  skin 
temperature  can  be  found  using  the  thin  skin 
approximation.  The  derivation  leads  to  equation  (2.2) 


which  is: 


temperature  at  end  of  jth  time  step 

(K) 

temperature  at  beginning  of  jth  time 
step  (K) 

time  step,  1  second 

temperature  of  ambient  air  (K) 

Cppd  (J/m2-K) 

h  =  convective  heat  transfer  coeffient 
( J/rn  -s-K) 

2 

Q  =  absorbed  thermal  energy  (J/m  ) 

The  missile's  velocity  and  altitude  are  extremely 

high  as  it  comes  within  range  of  the  space  based  laser 

(Bailey,  1985).  Under  these  conditions,  the  local 

convective  heat  transfer  coefficent,  h,  will  be  very 

small.  Using  a  procedure  explained  by  Hall,  h,  can  be 

calculated  (Hall,  1985:85-87).  For  example,  at  100  km 

altitude,  and  3000  m/s  velocity,  h  is  approximately  .04 

J/m  -s-K.  Then  using  equation  (C.l),  and  substituting 

in  the  appropriate  values,  the  radiation  absorption 

6  2 

term,  AQ.  is  approximately  1.5  *  10  J/m  . 

This  value  was  checked  by  two  seperate  methods. 

The  first  was  taken  from  notes  on  laser  weapons  by 
Bailey  (Bailey,  1985:45-46).  This  analysis  determined 
the  amount  of  energy  needed  to  melt,  or  remove  a  depth 
of  surface  in  a  given  amount  of  time.  Neglecting  loss 
process,  and  assuming  all  the  energy  goes  into  melting 
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where 


AT  = 

T  .  = 

air 

a  = 


the  material,  the  absorption  terra  is  approximately  3.0 
6  2 

10  J/m  .  The  second  method  was  taken  from  Raedy,  and 

involves  a  two  dimensional  heat  flow  equation  to 

calculate  the  temperature  at  any  point  and  time  in  the 

laser  spot  (Raedy,  1971:75-85).  This  analysis  also 

involved  melting  the  material,  and  the  calculated 

6  2 

absorption  term  was  approximately  2.75  *  10a  J/m  . 

Since  complete  melt  through  is  not  required  to  insure  a 
kill,  these  two  estimates  are  considered  high.  However 
some  melting  surely  occurs  prior  to  failure,  so  for  thi 
study,  the  minimum  absorbed  energy  density  was  assumed 

t-o  ho  9  fl  ♦  1  T/m 


2 


Appendix  D.  Data  for  Figures  i n  Chapter  IV 

This  appendix  contains  tables  of  data  used  to  plot 
figures  7,  8,  12  and  14. 


TABLE  D-I 

DATA  USED  TO  PLOT  FIGURE  7,  MODEL  COMPARISON 
OF  8,  16,  AND  32  MODELS 
PROBABILITY  OF  SURVIVAL  VS.  ROTATION  RATE 


Omega 
( r/s ) 

Model 

8 

C  nodes ) 
16 

32 

0.0 

.495 

.495 

.495 

0.1 

.728 

.635 

.618 

0.2 

.897 

.897 

.868 

0.3 

.994 

.993 

.991 

0.4 

.999 

.999 

.999 

0 . 5 

.999 

.999 

.999 

Missile  #41 

Launch  time:  0  sec 

Time  of  first  burst:  0  sec 
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TABLE  D-I I 

DATA  USED  TO  PLOT  FIGURE  8 


Omega 
( r/s ) 


MaxT 

Omega 

MaxT 

(K) 

(r/s) 

(K) 

708 

0.9 

467 

691 

1.0 

477 

651 

1.1 

472 

603 

1.2 

473 

555 

1.3 

474 

513 

1.4 

465 

487 

1.5 

452 

465 

1.6 

446 

457 

Missile  #41 
Launch  time:  0  sec 
Time  of  first  burst:  0  sec 
Fireball  rise  considered 


TABLE  D-I I I 
DATA  FOR  FIGURE  12 
4-0N-1  RANDOM  ATTACK 


Omega 

Ps 

Omega 

Ps 

(r/s) 

( r/s ) 

0.0 

.733 

0.9 

.259 

0.1 

.547 

1.0 

.475 

0.2 

.371 

1 .  1 

.637 

0.3 

.382 

1.2 

.725 

0.4 

.286 

1.3 

.829 

0.5 

.  144 

1.4 

.884 

0.6 

.080 

1.5 

.946 

0.7 

.090 

1.6 

.973 

0.8 

.143 

Missile 

#41 

Launch 

time: 

0  sec 

Time  of 

first 

burst:  0  sec 

TABLE  D- 1  V 

DATA  FOR  FIGURE  14,  SPOT  SIZE  COMPARISON 
2000  J/CM*2.  10,  20,  30  AND  40  CM  RADII 


Iraega 

Spot  Size 

r /s ) 

(cm) 

10 

20 

30 

.0 

0.00 

0.00 

0.00 

.  1 

0.00 

0.00 

0.00 

.2 

0.00 

0.00 

0.00 

.3 

0.07 

0.00 

0.00 

.4 

0.47 

0.00 

0.00 

.5 

0.98 

0.01 

0.00 

.6 

0.99 

0.07 

0.00 

.7 

0.99 

0.47 

0.01 

.8 

1.00 

0.98 

0.03 

.9 

1.00 

0.99 

0.07 

.0 

1.00 

0.99 

0.07 

.1 

1.00 

1.00 

0.49 

.2 

1.00 

1.00 

0.53 

.3 

1.00 

1.00 

0.98 

.4 

1.00 

1.00 

0.99 

.5 

1.00 

1.00 

1.00 

.6 

1.00 

1.00 

1.00 

Appendix  E .  Computer  Program 


This  program  is  a  modified  version  of  a  program 
called  "Therm",  written  as  part  of  a  master's  thesis  by 
Lt.  Barbara  A.  Hall.  The  program.  Therm,  was  written  in 
Fortran  77,  and  required  a  mainframe  computer  for 
effective  operation.  However,  several  simplifying 
conclusions  were  made  by  Lt.  Hall,  and  incorporated  in 
this  program.  Consequently,  the  following  program  is 
written  in  BASICA,  an  IBM,  and  Z-100  compatible  Basic 
language.  The  program  was  run  on  Z-150  personal 
computer. 


>  v  ^  r  v  V  V  t  v  V  :r  v  r  V  V  V  ¥  r  V  V  ¥  t  t  t  v  •)!  r  ^  n*  ^  t  V  t  ^  v  t  V  ^  i>'  v  v  #  *r  ^  t  v  •1'  n*  'i1  v  <i*  t  * 


3U"  -  *  *  *  * 

40  ’***- 


-k  3 

50 
60 
70 
30 
90 
100 
no 
120 
130 
140 
150 
160 
170 
130 
1  90 
200 
210 
-i  -> ,-, 

230 

240 

250 


*  *  vs* 

*  «  *  *  *  *  : 

f 


360 
370 
380 
390 
400 
410 
420 
430 
440 
4  50 
460 
470 
480 
490 
500 
510 

5>?r 

50Y- 

540 


Cumu 

lati 

ve  Ef  fee 

t  -  W  i  t 

h  Rotation 

Test 

Ver 

sion  #2 , 

16  nod 

e,  T-TMAX/2 

a*  a.  a>  a. 
t  -r* 

Prog 

-j-  a.  a*  u. 
t*  *r  -r  v 

ram  call 

a..<.a.^a.a.a.a> 

ed  1 T 1 6 

*<•  -a*  .la  »v  a.  a>  a. 

TT  T  T  T  T  T 

2  * 

a*  «u  «i»  u*  *i.  a<  <■<  ««r  .i»  u.  a. 

* *  £  * 

$$*** 

a.  ab  .a  a. 

a>  a.  a.  a« 

vU  u#  a#  a*  a#  a*  a#  »u  a«  «a  a#  a#  a#  a/  a*  a*  a*  a#  a*  a#  a#  %l  ^  «i« 

DECLARE  PROGRAM  CONSTANTS 

£  S  £  £  * 

a*  a.  a-  a-  a. 
'i*  'a*  *t*  t6  't* 

PI-3 

.141 

59: 

REM 

Radians 

R=1 1 

6.84 

“ 

REM 

Radius  of 

missi 

le  ( c 

THTA 

=0: 

REM 

Initial  an 

gular 

d  i  sp 

H0B= 

0: 

REM 

Height  of 

burst 

( km) 

ISS- 

619: 

REM 

Intensity 

Sure 

-Safe 

1 3  K  = 

S09 : 

REM 

Intensity 

Sure 

-Kill 

CE?  = 

200: 

REM 

Circular  E 

rror 

P  r  o  b  a 

a.  a.  .<•  «. 


(Kelvin, 


(Kelvin) 


•  v  v  'i*  'i*  *r  v  v  v  *r  *•*  *•*  *>  ■  ^  v  v  v  v  v  v  •»*  *!*  «i»  v  v  •«*  V  v  t  *»•  V  *i»  *}»  V  »|S  •[»  V  *1*  *{•  v  V  5}l  5,1  <,»  »Jt  J|»  v  *,•  V  •?  »?  »}»  *t*  * 


260  ’ 

-  DECLARE 

MISSILE 

SKIN  CONSTANTS 

;$>: 

1-  V  r  T  -i1  T  T  'I'  V  V  1*  Y  T*  'I*  V  T  T  V  1*  T  T  r  T  T  *1* 

V  T  ^  ^  *P  T  1* 

290  • 

300 

ALFA= . 5 : 

REM 

Absorptivity 

310 

CV  =  900 : 

REM 

Specific  Hea 

320 

D-.001 : 

REM 

Skin  Tnickne 

330 

RH0=2700 : 

REM 

Density  of  A 

340  ’ 
350  ’ 

^  a#  %v  a.  j.  a.  a,  a.  a.  .i.  ^  a.  a.  ..  a.  a.  w  a>  .«  a>  <v  a. 

*#*#$** 

a.  a.  ^  a.  a.  a<  a.  a.  >u  .i.  a.  a.  .u  a.  a.  < 

Aluminum 
lapacity  (v 
(m) 


( kg/m**3 ) 


•  a<  x  •«»  *J«  a>  a-  . ■•  .u  a.  j.  .•.  a.  a>  >u  a.  u.  «l  u,  . 


.  .t.  j,  .<•  a.  .u  a.  J.  a.  a.  a.  a.  a,  a.  a.  x 


DIMENSI0NALIZ2  STORAGE  ARRAYS  FOR  PROGRAM  CALCULATION 

>  a.  a«  a<  4>  a<  a>  >a  «'<  a*  a.  a,  a«  •>•  «w  a.  a#  a>  a.  >i>  <v  a*  a.  j,  a«  .u  a#  j,  a«  a>  a>  j.  %l  a>afa.M<a>a<xxxa>a>a<  »u  a.  a.  «l  a«  jl  a#  a*  a#  .l  , 


DIM 

DIM 

DIM 


VDATA(52),ZDATA(52) tXDATA(52) , ANGDATA ( 72 ) ,SX( 100) ,SY( 100) 

DCT(6,32) ,HFB(6, 32) f VEL(32) , ALT(32) ,DRD(32) ,ANG(32) ,H(32) ,TEMP(32) 
T2( 16) ,D0(16) ,RF( 16) .THETA ( 16) ,T1( 16) 


»  ■ 1  -  aj  a#  *i«  a#  a#  •'*  a«  *<j  a#  «<•  a*  a#  a#  a#  a#  «r#  .v  a*  »'r  a«  j#  «'•  % 

•  r  'l'  1*  •!>  "r  T  T  ■!*  T  T  •!•  T  T  *P  T  T  T  T  T  'I*  T  ■!■  -f  T  v  T  T  ‘ 


!  *  *  *  *  *  *  *  S 


INITIALIZE  ARRAYS  WITH  MISSILE  FLIGHT  DATA 


•  t  •••  *i*  i*  *r  w 


•  a.  -j,  a.  a.  a.  **.  a.  a«  «u  *’•  a.  -t.  a.  «■. 

*  'T*  *r*  't'  -r-  "f  'r'  *T*  o'  'r*  'r'  >T*  T*  T*  'T*  -T'  'i-  -i'  -T- 


*  *»*  •<*  *t»  «t»  •" 


a>  a.  a.  a.  . 


GOSUB  4490 
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INPUT  ATTACK/ 


PARAMETERS 


FOR  MN-41  TO  41 
BLUE-0 

INPUT  "ENTER  MISSILE  NUMBER" ; MN 

INPUT  "ENTER  DESIGNATED  LAUNCH  TIME" ; LT 

INPUT  "ENTER  NUMBER  OF  WEAPON  BURSTS  AFFECTING  MISSILE" ;MAXB 
INPUT  "ENTER  TIME  OF  FIRST  BURST  (  >=  launch  tine  and  even)";T31 
INPUT  "ENTER  YIELD  (kt)";Y 

INPUT  "ENTER  1  IF  YOU  WANT  FIRE3ALL  RISE,  0  IF  NOT"; RISE 
INPUT  "ENTER  ROTATION'  RATE  OF  MISSILE  (in  radians )"; OMEGA 
GOSUB  5300 
OMEGA-BLUE/20 

%  4  ,*•  ,  V  i  ^  %l#  4  4  *'#  *1*  -  -  ,1,  O#  %lj  ,*j  J-  •'#  *U  >1,  •'«  ,1,  *1*  -  -  -  ■  l#  «J*  ,U  hU  U#  4*,  4*4  4*a 

-r*  -»*■  *V»  »r»  »r*  'P*r“r-r'r*'r*Ti*,r,i,'i,'T,'i*i''PTT'<,'r'rv-rT'r'T'TTT'T*TT'i*'r'i*'.''r'rT'i,'i>’'’”r'r,TV,TT,i'TTT'i*'r'i''PT*i'i'T',' 

4I4  4  *  J4  J#  4*4  4/4  4^4  4I4  4*4  4I4  4&4  4i4|4  4^4  '  *4  4*4  4*4  4*4  4^4  444  4r4  4I4  4 14  41 4  44  4I4  4I4  4<4  4J4  4*4  4*4  4^4  4(4  4I4  4^*  4&4  4L4  4(4  4^4  4I4  4^4  4^4  4*4  4X4  4(4  4X4  4J4  4X4  4*4  4^4  4X4  4X4  4X4  4X4  4*4  4X »  4X4  4X4  4I4  ^4  ~  *  ~  4X4  4/4  4^4  4X4  -  X  -  ^ 

'f'l’  T'Tr,1**l''l**(4-(«TT,l''l,'l>*i>'I'TT'rT*r,T>1*T1''l‘T'r'r'*l*TTTTT'l'T'r**I''l*-T*'r'VTT'T*'l,T'^T'r,T'T'T'l’T1*TT'rTT,T'^T 

-  PRELIMINARY  PROBABILITY  DENSITY  FUNCTION  CALCULATIONS 
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3ETA-C1 ! /4 . 1 08 ) *L0G (ISK/ISS) 
ALPHA-. 5*LOG( IS5*ISK) 


PRELIMINARY  MATH  CALCULATIONS  FOR  THIS  ITERATION 


A«CV#RH0*D 

:  REM 

(j/m**:-;c> 

r  f=  .  is 

:  REM 

Thermal  fraction 
(surface  burst) 

XM-5 . 5 

:  REM 

Critical  point  along 
missile  length  (m) 

TMAX-.0417-Y* .44 

:  R  E i 

Thermal  Maximum 

BEGIN  CALCULATIONS  FOR  THE  MAXIMUM  TEMPERATURE 


CALCULATE  INITIAL  MISSILE  CHARACTERISTICS 


\T„<  1070  G0SU3  2340 

1030  ’ 

N  1090  '  ***>C:**********#iJt*#*:jE**:St:*:jES)<*:*sjs*******:**3*S:t:*#:*3t$£:***j*afs*#$:*:it****$*:)t****!|(>lt 


& 


1 1 ! 


W. 


1 1 

1130 
1140 
1150 
1  loO 
1170 
1  ISO 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
13  20 
1330 
13  40 
1350 
1360 
13 
13 


1390 
1400 
14  10 
1420 
1430 
1440 
1450 
1  4o0 

1  /  'J 

1  4S0 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
162$. 
1 6  'S.'-. 
1640 
16  50 
1660 
16  70 


CALCULATE  FIREBALL  CHARACTERISTICS 


U-Y/1000 
GOSUB  2620 


I  V  *  *  *  $  *  : 

1 


INITIALIZE  MAX  AND  MIN  TEMPERATURES  FOR  CALCULATION'S 

*  *p  *i*  t*  *i"  t  ^  t  *(•  'rV*r4*'r')*'r'r*r,*rT*,rVi,,r  *r  t  V  *v*  'i'  *}!  *{*  • 

LOV;T=  1 000000  !  :  H I T = 0  !  :  SUMPD=0! 


•  n*  «*»  «r*  »i®  *r*  'i'  *»•  f  *r*  * 


BEGIN  ITERATION  PROCESS  FOR  MAX  TEMPERATUR! 


*/*  t*  t*  t-  - 


*  *i*  *r“  *r  *1*  *i»  * 


GOSUB  1630 


'  T  *1«  'r  'I*  T-  T»  »■('  *T»  «(*  '|»  V  *■>»  *i»  ^  *1*  * 
»  *■*  *>'#  . '»  *1#  *U  •’*  Op  fc1 .  %*■  «  ** »  4*  »'*  « 

»  *»«  n*  <r  -n  **r»  •v  *i*  * ( »  ».»  *»»  -o  **•  *v  -c  *,»  • 


PRINT  THE  RESULTS 


'  -r*  '/*  *i* 


»  -r*  *r*  "i*  *•*  *r 


-r  *r»  *r»  *t»  f*  *\- 


LPR I 

NT : 

;  PRINT 

L?RI 

NT 

"MAX  TEMPERATURE 

(kelvin) 

"  ;  HA 

XT 

LPR  I 

NT 

"MISSILE  NUMBER" 

;MN 

f  D3  T 

o  t  u  1 

.i  i 

"ROTATION  RATE  ( 

radians)" 

;  0ME 

GA 

LPR  I 

N'T 

"PROBABILITY  OF 

SURVIVAL 

(Ps) 

»  . 

» 

-U  O-  .Ip  .4. 


i  i  i  i  *^  *  *'•  ^  *l»  *  J*  »•»  »■*  *'•  >i«  *i#  »>»  wp  *i.  * 


•  -i*  »*»  *r*  “r  t»  -r*  - 


*  *  £  *  S 


*  u.  *  j.  x  *  «u  *  ^  a  a  .i.  ...  .i,  .i.  4.  ...  *  .■,  *  *  *  *  *  .i.  o,  4.  *  j.  .y  .i,  »u  *’»  **»  o.  *v  op 

T'l'T'r.'fTT'i’  ‘r,rT1"TtTTTTT*rfTTTT^'rYTrT<M*TV1*'I*'('T*n*1' 

TRY  ANOTHER  MISSILE/BURST  SCENARIO  OR  QUIT 

P  »'j  »'p  Op  J»  i'p  .1/  «U  kC  p’p  U;  *U  *'p  pip  v P  » 'p  Op  Up  pip  Up  Up  Up  Up  X  X  X  V#  Op  Op  Op  Op  X  Op  Op  pL  X  p’p  Op  Op  X  Op  Op  X  Op  .  .p  Op 


3LUE  =  3LUE+  1 

’  INPUT  "WOULD  YOU  LIRE  TO  TRY  ANOTHER  SCENARIO  (YES  OR 

IF  BLUE  <=  32  THEN  GOTO  700  ELSE  GOTO  1590 


NO) ?" ;TRY$ 


IF  TRY$="YES"  THEN  GOTO  600  ELSE  END 


.<  i 

END 

*  i 


Op  Op  Op  Op  X  .Ip  v’p  Op  O.  Op  Op  Op  O.  -■  •  Op  Op  Op  Op  Op  Op  .  .  O.  Op  Op  •'.  ...  X  Op  >1,  O,  .1.  ..  ,«»  O.  Op  .1.  Op  Op  Op  .1. 
******  *«*  *f*  ***  *r  **»*»*•*»•  t  pi*-  «x  *t*  *v»  *v»  ***  **»  *o  *r*  •»»  •**  *i»  *r»  *r*  ***  *r»  **»  *»*  ».•  »x  *.*  »i*  v  -i'  t  t  i*  r  -r  **. 


Op  X  0»  Op  Op  *^p  O.  Op  Op  Op  Op  Op  Op  . 


SUBROUTINE:  TCALC 

p  Op  Op  Op  O.  ...  Op  ...  O.  Op  Op  Op  Op  .’p  X  .y  X  X  v.  .1,  X  »'p  Op  Op  X  X  «'»  .Ip  O-  Op  Op 
*  *>*  *T*  “f*  T  T  -I*  *.*  »T  'l%a'  'P  *,*•  PT  'I'  *r>  V  *r  ’I*  T  *P  'I*  •»*  P|»  'I*  P|"  p».  -y. 


To’ 


/  V'  •.'  •."  V' 


» 


EHS 


! 


«v> 


1700 
1710 
1720 
17  30 
1740 
1742 
1750 
1760 
1770 
1750 
1790 
1800 
1810 
1320 
1530 
1540 
15  50 
I860 
1370 
1380 
1390 
1900 
1910 
If**) 
1*60 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 

/*) 

2050 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2^ 
2180 
2190 
22 


N  B-l 
J®0 

T1=TEM?(J) 

T2=T1+1 ! 

I 

FOR  L=0  TO  15:  T1 ( L ) =TEMP( J ) :  NEXT  L 

I 

IF  ( T2  >  =  T 1  Ai.'D  J  <=  20)  THEN  GOTO  1770  ELSE  GOTO  2170 
J=J+1 
COUNTER® J 

IF  (NB  <  HAX3  AND  J-TMAX/2  >=  2*N3)  THEN  GOTO  1300  ELSE  GOTO  1330 
N3=NB+ 1 
NEWEST® 1 

GOTO  1860 

N  EWBST=0 
'  END  IF 

f 

FOR  L=0  TO  15:  DQ(L)=0:  NEXT  L 

t 

FOR  K-l  TO  N3 
TB=TBl+2*( X-l ) 

S3=(TB+2 ) / 2 

1 

f 

IF  ( NEWEST® 1  AMD  K-NB)  THEN  GOTO  1940  ELSE  GOTO  1990 
T-TBl-LT+J*TMAX/2-(TBl+J*TMAX/2-T3)/2 
GOSUB  3620 

VEL=V:ALT=Z:DRD=X:ANG=PHI 
GOSUB  3510 

GOTO  2020 
'  END  IF 

I 

GOSUB  3960 

? 

FOR  L=0  TO  15 

DQ(L)-DQ(L)+DC7(X, J)*TF*Y*TAU*CF*RF(L)*4. 1 36E+1 2/ ( 4 ! *PI* ( SR  *  2 ) ) 


T2(L)=(Tl(L)*(A-H(J)*THAX/4)+H(J)*THAX/2*TEH?(J)+ALF A*DQ(L))/(A+H( J)*TH 


NEXT  L 
GOSUB  5120 
NEXT  K 

FOR  L=0  TO  15 

IF  T2(L)  >  T 1 ( L)  THEN  T1(L)=T2(L) 
NEXT  L 
T1=T2 
GOTO  1760 
HAXT-Tl 

IF  ( MAXT  <  LOUT)  THEN  LOWT-MAXT 
IF  (MAXT  >  HIT)  THEN  HIT-MAXT 


4 


r  r 


22  1  Q 

Z  =  (  L0G(;iAXT )  - 

ALPHA 

2  2-*N 

ZP=ABS ( Z ) 

2230 

PZ=1 ! - 1 ! / ( 2 ! * 

(  1  !  +  • 

2240 

f 

2250 

IF  (Z  >=  0! ) 

THEN  ; 

2260 

r 

2270 

t 

2280 

r 

2290 

SU.MPD-SUMPD+PDI 

2300 

1 

2310 

1 

2320 

RETURN 

2330 

f  J.  .1.  J.  d.  J>  a- 

•A.  ^  -A.  Or  . 

T  T  T  1*  1 

2340 

’  **:******** 

«L.  •<.  ~U  -4-  vL  . 

T*  ***  *»“'»*  ' 

2350 

•  * 

2360 

»  £ 

2370 

1  * 

2380 

***** 

2390 

J  =  0 

2400 

T0=TB 1 -LT 

2410 

T=T0 

2420 

GOSUB  3620 

2430 

VEL( J)=V: 

ALT( 

2440 

1 

2450 

f 

2460 

GOSUB  4160 

W 

2. 

H(J)-H;  TE.’I?(  J 

» 

2490 

'Calculate  mi 

ssi  le 

2500 

'of  each  time 

step 

2  510 

» 

2520 

FOR  J=1  TO  31 

2530 

TM«T0+( J- 

.5)*T 

2540 

t-t:i 

2550 

GOSUB 

3620 

2  560 

VEL  ( 

J )  =  V : 

2570 

GOSUB 

4160 

2530 

H  ( J ) 

=  H : 

2590 

» 

2  600 

NEXT  J 

2610 

RETURN 

2620 

1  *  «V  4>  «)..'>  «>  «V  *'■*  •Jm 

it  **»  -r*  ~c  t 

2630 

'  * 

:p*4) 


•V‘r^*'r*r*i**rv¥'^''*VT'vv*r*PV*f»*!»» 


SUBROUTINE:  INITCHAR 


.  .1.  .1,  j,  »•«.  -•-  -J,  .1.  .1,  u-  %>«  ><•  *•»  o> 

*  -T"  *r*  'V  »n  't*  n*  *r*  *t»  *.*  *v*  t*  *<»  **»  «•«» 


SUBROUTINE :  FIRE3ALLCALC 


2c,<n 


2  o  o’u 
2670 
2680 
2690 
2700 
2710 
2720 
2720 
2  740 
2  7  50 
2760 
2770 
2  7S0 
2790 
2800 
2310 
2820 
2330 
2840 
2850 
2860 
2870 
2380 
2390 


FOR  \ = 1  TO  4 

FOR  J=1  TO  11 
DOT ( X , J ) =0 ! 
HF3(K,J)«0! 
NEXT  J 
NEXT  K 
J  =  0 

s-o 

NB=  1 


•FOR  EACH  TINE  STEP  J,  FIND  DIFFERENTIAL  FLUENCE  AND  FIREBALL  HEIGHT 


FOR  J-l  TO  31 

IF  (No  <  i  1 A X B  AX'D  J*T:iAX/2  >=  2*N3)  THEN  GOTO  2320  ELSE  GOTO  23 
N3=N3+ 1 
NE’.v  3ST  =  1 
GOTO  2360 

NE'.:3ST=0 


’FOR  EACH  BURST  X  THAT  HAS  OCCURED,  FIND  THE  THERMAL  FLUENCE 


2900 

29ns 

29Y0 


2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3lJjs* 
3  lVj 


FOR  X-l  TO  N 3 

TB-T31+2*(K-1) 

T  ?=  ( T  3 1  +  J  =>-  TMAX  /  2  -T3 )  /  TM  A  X 

IF  (TP  <=  10)  THEN  GOTO  2950  ELSE  GOTO  3070 
G0SU3  3240:  CTU-CT 

TP=TP- . 5 


IF  (T?  >  0)  THEN  GOTO  2990  ELSE  GOTO  3020 
G0SU3  3240:  CTD=CT 


GOTO  3050 
CTD=0 ! 


DCT(X, JJ-CTU-CTD 


Find  fireball  rise  for  burst  k  at  tine  t 


IF  (RISE  =  1)  THEN  GOTO  3110  ELSE  GOTO  3200 

IF  (NEW3ST-1  AND  X-N3)  THEN  GOTO  3120  ELSE  GOTO  3140 

t»(tbi+j*t;:ax/2-tb)/2! 

GOTO  3170 

T*TBl+( J-.5)*TMAX/2-TB 


HFB(K, J)-21 640.8* (W*. 177) *( 1 ! - ( 1 ! -T/240 ! )A 2 ) 


SBsSB&JSR 


3  ;v.‘j 
*■  b  » 

3  1  VO 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3230 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3350 

3390 

3400 

3410 

3420 

3430 

3«f0 

34o0 

3460 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 


.1  l  X  T  SC 
NEXT  J 

f 

2ETURX 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  *1*  V  ^  V  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  *1*  V  *!»  «{•  *i»  3}!  *{»  »»»  S{S  5}»  S»»  5{S  J}i  SjS  5|C  Jy!  Jjs  S{*  *!*  Si* 

t 

'  *  SUBROUTINE :  TliERMCOEF  * 

*  * 

» 

IF  (TP  <=  .75)  THEN  CT— . 02*TP+ . 24*( T? *  2 ) 

f 

IF  (T?  >  .75  AND  T?  <=  1.5)  THEN  CT= . 32*TP- .  1 2 

f 

IF  (TP  >  1.5  AND  T?  <=  2.5)  THEN  CT—  .  2 5 7 2 1 9+ .  5564 1 5*TP- .  0969029- ( TP A 2  ) 

» 

IF  (TP  >  2.5  AND  T?  <  10)  THEN  CT= . 335808+ . 0949904*T?-4 . 945 1 4E-03* ( TP  *  2 ) 

f 

IF  T ? >  = 1 0  THEN  CT= . 8 


RETURN 

T  **.  >•<  %<u  .*»  ai.  w.  ai.  .1.  ~u  j.  %u  <j< 

•i»  *i*  t  *»•  *r  *r  'I*  *r  v  *r  *r  *r  'f  *t* 


a*.  »l»  J.  J.  aU  J«  4.  aW  J.  »•*  *•»  vU  aV  a*.  aV  **•  **»  «**  *t*  «*»  J.  J»  *1*  «U  aU  .u  .U  %'•  .<•  at.  «J«  4.  J.  J.  a».  J.  « 


SUBROUTINE:  ALTMISSPOSN 

Used  if  burst  occurs  during  time  step,  and  new  missile 
position  characteristics  required.  If  not,  gosub  5760. 


•  T  T  T  1"  V  T  'I*  T  -r  T  1"  T  V 


DELTAZ 
X  3  =  S  X  ( 
GR  = 
SR» 


TAU 

GOSUB 

RETURN 


=ALT-HFB(X, J) 

S3):  Y3=SY ( SB ) 

3QB ( (  SY  (MIJ )-YB )  *2+( DRD-X3)  * 2 ) 

SQR ( GR *  2  +  DELTAZ  *  2 ) :  SRX=DRD-X3:  SRZ=DELTAZ 
IF  VEL=0  THEN  GOTO  3560  ELSE  GOTO  3570 
CF= 1 :  GOTO  3590 

COSPHI=(SRX-VEL*COS(PHI)+SRZ*VEL*SIN(PHI) )/(SR*VEL) 

CF=SQR( 1 ! -COS PHI *2) 

=EX?(-.02455-o.439E-05*SR-1.407E-09*SRA2+l . 792E-14*SR*3 ) 

4560 : COUNTER=COUNTER+ 1 

a/»  w.  a1.  J.  O.  .<>  J;  V.  «*»  .*>  J,  aL  alf  Ju  U.  al.  a^  afa  4.  a'<  ab  ab  ab  aU  al,  a',  ab  >V  «■>  »'»  a1,  ab  a>.  a>,  J,  ai#  J.  >1.  J.  X  aU  a1.  a'«  a*,  af.  X  .1#  XXX  a',  X  a'.  J>  a',  ab  X  V#  X  X  X  X 


3630 
3Cy-2 
3  6\’~J 
3660 
3670 
3o30 
3690 
3700 
3710 
3720 
3730 
3740 
3750 
3760 
3770 
3780 
3790 
3800 
3810 
3820 
3830 
3840 
36  50 
5860 
33  70 
3330 
3890 


3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4030 

4090 

4100 

4110 

4120 

4130 

4140 

4UQ 

4 

4170 


*  SUBR0U7I .'.’3 :  .IISSILECALC 

T  V  V  V  V  T  T  V  r  T  t  ’l*  ¥  *!•  T  T  T  T  T  V  '•  If  j|:  ^  ^  ^  ^  ^  ^  5]:  ^  ^  ^  ^  ^  ^  ^ 


•  T  4  V  • 


TD-INT(T):  TU-TD+1:  REM  UPPER  AND  LOWER  INTEGER  VALUE  OF  T 

IF  T>50  THE:.’  GOTO  3710  ELSE  GOTO  3730 

V=105 ! *T-1350 
Z=2460!*T-65000! 

X=3700 ! *T-1 2 7000 ! 

GOTO  3320 

V-(T-TD)*(VDATA(TU+1)-VDATA(TD+1))+VDATA(TD+1 ) 
Z=(T-TD)*(ZDATA(TU+1 )-ZDATA(TD+l ) )+ZDATA(TD+l ) 
X-(T-TD)*(XDATA(TU+1)-XDATA(TD+1))+XDATA(TD+1) 


I?  (T  >  70)  THEM  GOTO  3340  ELSE  GOTO  3360 
PHI  =  30.6-.ls*:(T-70) 

GOTO  3330 

?iii-(t-td)*(a;;gdata(tu+i)-angdata(td+i  ))+a:;gdata(td+d 


rstur:: 

r 

»  *** 


V=V* . 3048 
Z=Z* . 304S 
x=x*.3048+sx(:i:o 
PHI-PKI*PI/180 


%•»  j.  j.  x  j.  4.  «>»  •>. 
*«»  *t*  *i«  t  *r  t 


* 


**** 


** 


***** 


u.  «« 
ir  *r*  -»* 


1  ±. 


’  *  SUBROUTINE:  MISSPOSN 

»  * 


■J*  ^  O.  X. 


» 

» 


*■»  X  V.  *•.  X 
"**  *T>  *T*  -T*  -T*  -r 


X  X  X  X  X  X  X  X  X  4, 


x  x  x  x  x  x  4. 
¥  «r  t*  -r  t  *(•  *r 


DELTAZ=ALT(J)-HFB(K,J) 

XB-SX(SB):  YB=S Y( SB  ) 

GR=SQR( ( SY(MN)-YB) *2+(DRD(J)-XB)*2) 
SR=SQR(GR*2+DELTAZ*2) :  SRX=DRD( J ) -XB :  SRZ=DELTAZ 


IF  VEL( J )=0  THEN  GOTO  4060  ELSE  GOTO  4070 
CF=1 :  GOTO  4090 

COSPHI=(SRX*VEL(J)*COS(ANG(J)  )+SRZ*VEL( J)*SIN( ANG(  J)  ))/(SR*VEL(  J) ) 
CF=S0R(1!-C0SPHI“2) 

TAU-EXP( 02455-6. 439E-05*SR- 1 .407E-09*SR* 2+1 .792E-14*SR *3) 
DP.D-DRD(J):  GOSUB  4860 : COUNTER-COUNTER* 1 


RETURN 

•  X  X  X  X  X  4.  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X 


*  #  *  £  *  £ 


*  *  £ 

*  * 

»  * 


*  %  *  sjt  jjt  *  * 


SUBROUTIN'E:  HEATCALC  (heat  transfer  coefficient  h  (j/m2-s-k) 


4  i.  ;„-n 

4  r'-'-o 
4200 
4210 
4220 
4230 


4250 


4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4 

4fti0 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4  5  70 
4580 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
0 
0 
0 
J 


,u  .k  ^  u.  .y  ■'»  •»'.  ..u  •>.  .•» 


Z=ALT ( J )  :  G0SU3  4310:  RE.’l  Find  ambient  air  parameters 
CP=240*4.184 

RE-RHOA*VEL( J)*XM/MU 
PR«iiU*CP/SA 

IF  RZ<=500000  !  THE:,'  NU=  . 332-PR “  .  333*RE *  .  5  ELSE  HU=  .  0296-P2  *  .  3 
H»:.:u*SA/x:-i 


SU3R0UTI 


Standard  Atmospheres  (  47  km  ) 


IF  Z<  1 1000  THE:,’  LX=-. 006545:  PK  =  101300  !  :  TX«283.15:ZK-0! 

IF  Z>=11000  AMD  Z<20000  THEN  LK=0!:  PK=22690 ! :  TK=216.65:  ZK=11000 
IF  Z>=20000  AND  Z<32000  THEN  LK-.001:  ?K=5528!:  TK=216.65:  ZK-20000 
IF  Z>=32000  AND  Z<47000!  THEN  LK=.0028:  PK=888.8:  TX=223.65:  ZK-32000 
IF  Z>47000!  THEN  PRINT"  Consult  NOAA  for  values" 

IF  LK  =  0  THEN  GOTO  4410  ELSE  GOTO  4420 

P-PK*EXP(-.034164*(Z-ZK)/TK):  TA-TK:  GOTO  4430 
TA=TK+LK*(Z-ZK) :  P-PK*(TK/TA) *( .034164/LX) ;  GOTO  4430 
RHOA».003484*P/TA 

M U  =  (  1 .458E-06*TA*1 .5)  / (TA+110.4):  REM  (lcg/m-s) 

:<A  =  (2.64638E-03*TA*1.5)/(TA  +  245.4*10a(-12!/TA)):  REM  ( j/m-s-k) 


RETURN 


SUBROUTINE:  INITIAL  DATA  INPUT 


^  "i*  t*  4"  V  t  *p  i*  *1*  n*  ^  V  v  *r  'J*  V  v  *J“  *r  V  v  v  v  v  v  t  "r  v  t  *?  m*  v 


OPEN  "I",r:l veldata.txt" 

FOR  1  =  1  TO  51:  INPUT  yl.VDATA(I):  NEXT  I  :  CLOSE  ,7 1 
OPEN  "I",  7l , "altdata.txt" 

FOR  1=1  TO  51:  INPUT  #1,ZDATA(I):  NEXT  I  :  CLOSE  #1 
OPEN  "I",? 3," drddata.txt" 

FOR  1=1  TO  51:  INPUT  ?3 ,XDATA( I ) :  NEXT  I  :  CLOSE  #3 
OPEN  "I",„ll ,  "degdata.txt" 

FOR  1  =  1  TO  71:  INPUT  ? 1  , ANGDATA( I ) :  NEXT  I:CLOSE  #1 
’OPEN  "I",#l ,"MLO.txt" 

'FOR  1=1  TO  100:  INPUT  #1,MLT(I):  NEXT  I  :  CLOSE  H 


«  «**«**««##««*  JSt:*:***:*#**:*****#:*:******:*#:*;*:*:,***:*!***#**#*;*: 

»  * 

*  *  SILO  POSITIONS  (X  &  Y  COORDINATES) 


*1#  sV  ^  J.  « 


4  710 


4  *"v': 

4  7"Jo 


4/4Li 

4750 

47o0 

4770 

4730 


Da =51 9. 62  :  DY=300 

FOR  1=0  TO  95  STEP  5 

S Y( 1+1 ) =0 :  SY(I+2)-2*DY:  SY(I+3)=4*DY:  SY(I+4)=DY:  SY(I+5) 

NEXT  I 


=  3*DY 


4310 
4S20 
43  30 
4340 
48  50 
4360 
43  70 
4880 
4890 
4900 
4910 
4920 


J  =  G 

FOE  1=0  TO  38  STEP  2 

5X( J+1)=I*DX:SX( J+2)=I*DX:SX( J+3)=I*DX:SX( J+4)=(I+1)*DX:5X( J+5) 
J  =  J  +  5 

NEXT  I 

retur:: 

f  -v  »v  u.  «v  *•»  -y  >'<  >u  »•«  •>  >’<  a  .y  .i.  y,  «i>  «i.  ^  o>  %y  »'<  j.  >u  >>•  >'•  >'<  »■-  w  .'j  •'<  *•»  u#  «•«  .  -  *»,  «j*  *>, 


=(I+1)*DX 


•  j.  «v  «a»  »<u  •  •  «'•  -a.  j,  »•.  j.  -u  u.  .a.  «.i#  «•.  *j.  .x*  «y  «y  »y  «'•  «v  «u  »t,  »y  «v  •*»  a  .i.  »•<  >u  »■-  *y  >;<  a.  <■<  »l.  »y  ,v  v#  *y  ^  ^  %u  •>.  x  «.'»  •'»  *  « 


ROTATION  FACTOR  CALCULATIONS 


«*•  v.  .v-  «i»  »».  o»  a.  a*  v.  «w  a>  «u  a*  «u  «*»•*»  «*»  *v  «*>  *u  »•»  •**  «**  «*-  J-  »*»  »■»  *a.  .«»  «j-  a.  «i»  •>.  -t.  *i.  a.  .u  %•»  >•<  -y  -j.  u.  j-  *•*  .<«.  -j.  -j.  -u  -j.  -x,  ju 


4940 

4950 

4960 

4970 

4flr* 


J-COUNTER  THE:;  GOTO  4910  ELSE  GOTO  4960 
■TA(0)»TH7A+THETA(0) 

FOR  1=1  TO  15 

THETA(I)»THETA(I-1)+PI*.125 
NEXT  I 


BOTTO:i=DRD-X3 :  T0P=SY(MN)-Y3 

if  3otto::=o  the::  gamma »o  else  gamma-atsctop/bottom) 


5000 

5010 

5020 

5022 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5100 

5110 


FOR  1=0  TO  15 


IF  DRD(J)  >  XB  THEN  THETAPRIME-( 1 . 5*PI )+GAMMA 
IF  DRD(J)  <=  XB  THEN  THETAPRI!1E=(  .  5*PI ) +GAMHA 
ZPRi:iE=ABS(THETA?RIME-THETA(  I ) ) 
RF(I)=COS(ZPRIME) 

IF  RF( I ) < . 25  THEN  RF(I)=.25 
THTA=0MEGA*TMAX/2 


NEXT  I 


;turn 


«*.  .v  j.  »■.  a  j»  u,  u,  »b  a.  «l  4.  «u  >■<  .i.  «».  »>.  »>.  »■»  *V  «■»  »v  9«  A  »!?  »!«  V*  •>  t!*  ;!»  j;  *!*  j*  A  ^  A  A  *  9?  9»  9*  « 


f 


PRINT  THE  RESULTS 


5120 

s  •r- ' 
*  *  *  . 

5  iVu 

5150 

5160 

5170 

5  1  SO 

5190 

5200 

5210 

5220 

5230 

5240 

5250 

5260 

5270 

5280 

5290 

5300 

5310 

5320 

5330 

5340 

5350 

5360 


T»  -*•  f  -T*  -t* 


i  V  T  T  V  V  T  V  T  ^  T  T  Y  T  V  V  V  ^  ^  V  T  V  T  V  V  V  >T  V  T  V  T  V  T  T  V  T  T  T  1*  T  ^  T  V  -I*  V  V  T  V  T  >T  T  T  “f  T  ¥  V 


' LPRINT  "TIME  STEP" ;J, 
'  LPR I  NT  "NODE  TEMP 
FOR  1=0  TO  15 


"BURST  NUM3EE " ; K , "T/TMAX" ; ( ( J- . 5 ) *TMAX/2 ) /TMAX 

THETA  ROTATE-FACTOR  GAMMA" 


IF  T2  (  I )  >  T 2 
' LPRINT  USING 
' L  ?  R I N T  USING 
'LPRINT  USING 
' LPRINT  USING 
'LPRINT  USING 
NEXT  I 

f 

RETURN 

t 

LT  =  0 
TB1=0:  MAXB=4 

Y=2000 :  RISE=1 
THTA-0 

FOR  1=0  TO  15:  THETA(I)=0:  NEXT  I 
RETURN 


THEN  T2=T2 ( I ) 

#«?»";  1  +  1, 

;T2(I)f 

4 i 4 . 4 4?"  ; THETA ( I )  *  180 /PI , 
44 .4*44" ;2F(I) , 

4 4 4 4 . 4 4 4 " ; G A MM A* 1 8  0 / P I 
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The  study  investigated  several  different  scenarios.  The  first 
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a  4-on-l  attack  of  a  missile  launching  system.  More  specifically, 
one  missile  was  subjected  to  four  bursts  located  in  various 
positions  surrounding  the  missile.  The  intent  was  to  determine  if 
rotating  a  missile,  even  when  surrounded  by  thermal  radiation,  would 
increase  the  probability  of  survival.  Finally,  the  missile  is 
attacked  by  a  space-based  laser  with  a  maximum  absolute  power  of  10 
megawatts.  In  all  cases,  the  rotation  rate  was  limited  to  a  maximum 
of  1.6  radians/second,  as  established  by  studies  at  the  Air  Force 
Institue  of  Technology.  Using  computer  programs,  the  maximum  skin 
temperature  was  calculated,  with  the  resultant  probability  of  damage 
determined  using  a  cumulative  log-normal  distribution  function. 
Comparisons  were  made  between  the  rotating  and  nonrotating  missiles 
to  determine  if  rotation  did  increase  the  probability  of  survival 
for  the  missile  system. 

In  all  scenarios  studied,  rotation  significantly  decreased  the 
maximum  skin  temperature,  increasing  the  probability  of  survival  for 
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amount  of  protection.  As  a  consequence,  even  at  these  relatively 
low  rotation  rates,  rotation  is  an  effective  defense  against  the 
thermal  threat. 


